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ABSTRACT 

We present measurements of the color and luminosity dependence of galaxy clustering at z 
1 in the DEEP2 Galaxy Redshift Survey. Using volume-limited subsamples in bins of both color 
and luminosity, we find that: 1) The clustering dependence is much stronger with color than with 
luminosity and is as strong with color at z ~ 1 as is found locally. We find no dependence of the 
clustering amplitude on color for galaxies on the red sequence, but a significant dependence on color 
for galaxies within the blue cloud. 2) For galaxies in the range L/L*~0.7 — 2, a stronger large- 
scale luminosity dependence is seen for all galaxies than for red and blue galaxies separately. The 
small-scale clustering amplitude depends significantly on luminosity for blue galaxies, with brighter 
samples having a stronger rise on scales rp < 0.5 Mpc. 3) Redder galaxies exhibit stronger small- 
scale redshift-space distortions ( "fingers of god" ) , and both red and blue populations show large-scale 
distortions in ^(rp,7r) due to coherent infall. 4) While the clustering length, tq, increases smoothly 
with galaxy color (in narrow bins), its power-law exponent, 7, exhibits a sharp jump from the blue 
cloud to the red sequence. The intermediate color 'green' galaxy population likely includes transitional 
galaxies moving from the blue cloud to the red sequence; on large scales green galaxies are as clustered 
as red galaxies but show infall kinematics and a small-scale correlation slope akin to the blue galaxy 
population. 5) We compare our results to a semi-analytic galaxy formation model applied to the 
Millenium Run simulation. Differences between the data and the model suggest that in the model 
star formation is shut down too efficiently in satellite galaxies. 

Subject headings: galaxies: high-redshift — galaxies: evolution — cosmology: large-scale structure of 
the universe 



1. INTRODUCTION 

The clustering of galaxies reflects both the spa- 
tial distribution of dark matter, which is dependent 
upon cosmological parameters, and the complex physics 
which governs the creation and evolution of galax- 
ies within their host dark matter halos. Cluster- 
ing measu res can therefore constrain both cosmol- 
ogy (e.g., iPeacock et all l2001t lAbazaiian et all l2005l : 
lEisenstein et al.ll2005f ) and galaxy formation models, in- 
cluding parameterized models like the halo occupation 
distribution (HOD) models (e.g., IBerlind fc Weinber 



distribution [tiUU) models (e.g., Ihierlmd &: Wembe] 
2005; 'Yan. Madgwick. fc Whiti Ipipol l^g et all [200 
Zehavi et al. 2005; P hleps et al.l l200a) or more phys- 
ical semi-an alytic models and hydrodynamic s i mula- 
tions (e.g.. ICen fc Ostrikeil l2000i : iPearce et all l2001l : 



[Weinberg et al.l 20041 ). In particular, the detailed depen- 
dence of clustering on galaxy properties such as lumi- 
nosity and color at different redshifts provides a wealth 
of information on galaxy evolution. The formation and 
evolution of dark matter halos is straightforward to pre- 
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diet and has been characterized with large N-body sim- 
ulatio ns (e.g.. lMo fc Whitdll996l : ISheth. Mo. fc TormenI 
l2001h . but the details of how baryons in galaxies pop- 
ulate their parent dark matter halos are less well un- 
derstood and should depend on the detailed physics of 
gas accretion, star formation, and feedback processes. 
Comparisons of the observed dependence of clustering 
on galaxy properties at different epochs with theoretical 
models and simulations can illuminate the relative im- 
portance of various physical processes involved in galaxy 
evolution. 

On scales of ~ 1—20 h^^ Mpc the galaxy two-point cor- 
relation function is roughly a power law, ^(r)= (r/ro)"'^ , 
and has long been known to depend on galaxy properties, 
with one of the strongest dependencies being galaxy type. 
The division of the overall galaxy population into two dis- 
tinct types has been seen clearly as a function of galaxy 
morphology, color, and spectral type, where galax- 
ies tend to be either early/bulge-dominated/red/non- 
star form ing or late/di s k-dom i nated/blue/star form- 
ing (e.g., [S trateva et aL 200T; 'Madgwi ck et all l2002t 
iBlanton et al. 2003a; Kauffmann et al. 20Q3|). Restframe 
color has been found to correlate well with both spectral 
type and morphology, with some scatter. In this paper 
we focus on the dependence of clustering on restframe 
color, as this quantity can be determined more robustly 
for large samples of galaxies at both low and high redshift 
than, for example, morphological type. 

The correlation between g alaxy type and cluste ring 
amplitude was first seen by iDavis fc Gelled ()1976f ) as 
a difference in the clustering slope of morphologically- 
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selected elliptical and spiral galaxies, with ellipticals hav- 
ing a steeper slope. This implied that ellipticals are 
more tightly clustered than spirals and that a neigh- 
boring galaxy of an elliptical was more likely be to an - 
othcr elliptical. This was confirmed by Dressier (1980') 
who found a tight relation between morphology and lo- 
cal galaxy density in clusters of galaxies, with ellipticals 
preferentially located near the centers of clusters. 

The color dependence of galaxy clustering has been 
measured more prec isely in subsequent s urveys of the 
local universe fe.g.. iLovedav et al.l Il995t iHermit et al.l 
119961 : IWillmer. da Costa, fc Pellegrinil Il998f l. with the 
most statistically-significant results provided by the 
latest large surveys, the 2dF Galaxy Redshift Sur- 
vey (2dFGR S) and the Sloan Digital Sky Sur- 
vey (SPSS) (INorberg etH\ I2002b|: IZehavi et al.l [200l 



Madgwick et al.l l2003al: iLi et al.l l2006al: IZehavi et al l 



idg- 

2005f) . Table 1 provides an overview of clustering results 



for color- and spectral type- selected samples from 2dF- 
GRS and SDSS as well as recent intermediate redshift 
surveys. 

In the 2dFGRS, galaxies are classified into early- and 
late-type usin g a spectral analysis iii ste ad of using rest- 
frame colors ()Madgwick et al.l I2002D . iMadgwick et al.l 
(|2003ar ) find that early-type galaxies are much more 
strongly clustered than late-type, with a relative bias 
within spheres of 8 Mpc radius of brei = 1-45 ± 0.14. 
INorberg et al.l (|2002b( ) investigate the luminosity depen- 
dence of clustering of early- and late-type galaxies in the 
2dF data and find that while there is a strong difference 
in the relative mix of early- and late-types as a func- 
tion of luminosity, at L > L* this mix is not driving the 
overall increase of clustering for bright galaxies. Both 
galaxy types display a clustering increase at higher lumi- 
nosities, with a constant relative bias between early- and 
late-types as a function of magnit ude at L > L*. 

T hese trends are con firmed bv iBudavari et al] ()2003l ) 
and lZehavi et al.l ()2004D using SDSS data. These authors 
find that both red and blue populations show a significant 
increase in the correlation length with luminosity, though 
within the blue cloud there is no significa nt change in 
the slo pe of ^(r) for the brighter galaxies. IZehavi et al.l 
()2004f l find that within the red population, however, the 
faintest galaxies (—19 < Mr < —18) have a steeper slope, 
which is likely due to a higher fraction of these galaxies 
being satellites of ce ntral galaxies in groups and clusters 
(jBerlind et al.l[2005h . 

Color dependence of galaxy clustering has been de- 
tected in a v ariety of intermediate redsh i ft (0.2 < z < 1) 
surveys (e.g., IC arlberg et al."1997', '2001'; I Shepherd et"aLl 
120011 : iFirth et al. 2002; Phleps et al. 200|). The most 
robust results to date at z ^ 1 are from the VIMOS- 
VLT Deep Survey fVVDS: ,Le Feyre et al.ll2005h an d the 
DEEP2 Galaxy Redshift Survev (|Davis et al.ll2003h . Us- 
ing VVDS data, drawing from a total s ample of 6500 
;alaxi es at < z < 1.5 over ~ 0.5 deg^, iMeneux et al.l 
measure the clustering of red and blue galax- 
ies, dividing the two classes using the observed color 
bimodality. They find that red galaxies have a higher 
correlation length and steeper slope than blue galax- 
ies at z = 0.8, and that the relative bias of red to 
blue galaxies does not depend significantly on luminosity 
(|Marinoni et al.ll2005f). The se results are consistent with 
those of ICoil et all (|2004bD . who use an early DEEP2 



data sample of 2200 galaxies at 0.7 < z < 1.35 in ~ 0.35 
deg^ to determine that galaxies redward of the median 
color of the sample, {B — R)o = 0.7, have a significantly 
higher correlation length and steeper slope than bluer 
galaxies. 

Unfortunately, previous studies at intermediate red- 
shift covered relatively small volumes and are dominated 
by large cosmic variance errors. Here we improve upon 
earlier z ^ 1 DEEP2 results using the completed DEEP2 
data set and volume-limited samples. We measure galaxy 
clustering in several restframe color bins, with a particu- 
lar emphasis on 'green valley' galaxies at the intersection 
between the red/early- and blue/late- type populations, 
and as a function of luminosity within the red and blue 
galaxy populations separately. 

The correlation between galaxy type and clustering 
strength is also reflected in measurements of the 'envi- 
ronment' or local gal axy overde nsity around objects of a 
given class. Locally, H ogg et all (j2003) find a strong rela- 
tion between local overdensity and galaxy color in SDSS 
data, on scales of both 1 Mpc and 8 h~-^ Mpc, with 
red galaxies favoring more overdense regions. They find 
little dependence of overdensity on luminosity for galax- 
ies in the blue cloud (though the brightest blue galax- 
ies do appear to be in more overdense regions) and a 
strong dependence on luminosity for red galaxies locally, 
with both the brightest and faintest red galaxies being 
in denser environments than intermediate luminosity red 
galaxies. The lack of a stronger luminosity-overdensity 
relation for blue cloud galaxies may be somewhat at odds 
with the clustering results from 2dF and SDSS, which 
find a steeper, roughly linear trend of relative bias with 
magnitude, in that brighter blue galaxies a re significantly 
more clustered than fainte r blue galaxies (INorberg et all 
I2002bl: IZehavi et al.l [20051) . ISwanson et al.l (|2007D . ho^ 
ever, also study the clustering of blue galaxies in SDSS 
and do not find a significant dependence on the rela- 
tive bias with luminosity, in the same magnitude range 
as these other local st u dies, in contrast to the results of 
INorberg et all (l2002br i: IZehavi et al.l (|2f)05h . 

Galaxy overdensity as a function of color and lumi- 
nosity was recently measured for DEEP 2 galaxies by 
ICooper et al.l (|2006( l. who find a similar dependence on 
color as is seen locally, and a dependence on luminosity 
for the brightest of the blue galaxies, which lie in over- 
dense regions. No dependence on luminosity was found 
for galaxies on the red sequence, in the luminosity range 
probed by DEEP2. We expand upon that study here by 
measuring the clustering as a function of scale (instead 
of using a scale-averaged overdensity) using a statistic 
that may more readily be compared with simulations, as 
it does not depend upon the bias of the tracer sample 
used, as is the case with overdensities. The two-point 
correlation function results presented here will be used 
for HOD modeling of the galaxy population at z ~ 1 as 
a function of galaxy type in a future paper. 

Significant deviations in the two-point correlation func- 
tion from a power-law form have recently been detected 
on scales r < 10 h^^ Mpc for both local and high 
redsh i ft populations (e.g IZehav i et al.' '2 0041 : iLee et all 
[20061: lOuchi et al.l [20M [CoireT al. 2006^ These de- 
partures from a power law are naturally explained in 
the HOD framework, which provides a statistical under- 
standing of the relation between galaxies and their dark 
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matter halos, as the transition between pairs of galaxies 
within the same halo (the 'one-halo' term) on small scales 
and galaxies in different halos ( the 'two-halo' term) on 
larger scales (jZehavi et al.ll20(>i) . Small-scale deviations 
from a power law are expected to be more significant 
at z > 1, where the exponential tail of the mass func- 
tion moves to smaller masses, so that at a given mass 
scale the slope o f the n iass function becomes steeper 
(iKraytsov et alj 120041: iCqnrov Wechsler. &: Kravtsovl 
l2006HZheng. Coil, fc Zehavill2007f) . Observationallv. de- 
viations are also seen at a given redshift for the brightest 
galaxy samples, with L > L*. This is likely due to the 
satellite fraction decreasing for brighter samples; brighter 
galaxies are mor e likely to be central galaxies than 
satell i te ga laxies (Zchavi ct al. 2005; Mandclbau m et al. 
2006t Ivari den Bosch et al. 2007; Zheng, Coil, fc Zehavi 



Details of the DEEP2 observations, catalog constructio n 
and data reduction can be found inlDavis et all (^2003^ . 



Coil et al.l 



Davis et al 



2007f ). which leads to a steeper one-halo term (Tinker 



et al, in prep.) as more pairs are central-satellite pairs 
rather than satellite-satellite pairs and have, on average, 
a smaller separation. 

Departures of the two-point correlation function for 
DEEP2 galaxies from a power law as a fu nction of lumi- 
nosity are presented in ICoil et al.l ()2006a[ ): here we per- 
form a similar analysis for galaxies in DEEP 2 as a func- 
tion of color. Given that the galaxy correlation function 
is not a perfect power law, quoted tq and 7 fits clearly 
do not fully capture all of the information present in 
clustering measures. However, they are still useful as an 
approximate measure of the amplitude and shape of the 
correlation function on scales of ^ 1 — 10 Mpc and 
facilitate simple comparisons between observational re- 
sults and models. As such, we provide ro and 7 fits here 
but also discuss deviations from a power-law form for the 
data samples investigated. 

The outline of the paper is as follows: §2 briefly de- 
scribes the DEEP2 survey and defines the data samples 
used here. In §3 we discuss the methods used in this pa- 
per to measure the two-point correlation function, infer 
the real-space correlation length and slope, and correct 
for observational biases. §4 presents our results on the 
color dependence of galaxy clustering at z ~ 1. In §5 
we compare our res ults with the semi-analytic models of 
iCroton et al.l ()2006f ) applied to the Millenium Run dark 
matter simulation, and we conclude in §6. 

2. DEEP2 GALAXY SAMPLES 

2.1. The DEEP2 Galaxy Redshift Survey 

The DEEP2 Galaxy Redshift Survey is a recently- 
completed projec t using the DEIMOS spectrograph 
(|Faber et al.ll2003l ) on the 10m Keck II telescope to sur- 
vey optically-selected galaxies at z ~ 1 in a comoving vol- 
ume of approximately 5x10^ Mpc^. Using ~ 1 hour 
exposure times, the survey has measured high-confidence 
(> 95%, quality flag Q = 3 or 4) redshifts for - 30,000 
galaxies in the redshift range 0.7 < z < 1.5 to a limiting 
magnitude of i?AB = 24.1. 

The survey covers three deg^ of the sky over four widely 
separated fields to limit the impact of cosmic variance. 
Each field is composed of 2-4 'pointings', each of which 
covers an area of ~ 0.5 by ^ 0.7 deg, corresponding 
to the area covered by the CFHT12k camera; there are 
10 pointings in total. The DEEP2 spectra have high 
resolution {R ^ 5,000), and rms redshift errors (deter- 
mined from repeated observations) are < 35 km s^^. 



2004al) iDavis. Gerke. fc NewmanI 1)2005^ and 
(|2007[ ). K-corrections, absolute Mb magni- 
tudes and rest frame {U — B) color s have been derived 
as described in IWillmer et al.l ()2006[ ) . We do not include 
luminosity evolution in the K-corrections for Mb- Abso- 
lute magnitudes given in this paper are in the AB system 
and are Mb — 5 \og{h) with /i = 1, which for the remain- 
der of the paper we simply denote as Mb- 

To convert measured redshifts to comoving distances 
along the line of sight, we assume a flat ACDM cos- 
mology with = 0.3 and r^A — 0.7. We define h = 
Hq/{1QQ km s^^ Mpc^^) and quote correlation lengths, 
ro, in comoving h^^ Mpc. 

2.2. Galaxy Sample Definitions 

From the full flux-limited DEEP2 data set we define a 
variety of volume-limited subsamples, each correspond- 
ing to some range in galaxy (U—B) color and Mb magni- 
tude. A volume-limited sample is not subject to redshift- 
dependent selection effects, i.e., a galaxy of given rest- 
frame properties will be within such a sample regard- 
less of its redshift. While volume- limited samples from a 
given survey will include fewer galaxies than flux-limited 
samples, they are much easier to interpret and compare 
with theoretical models, and are therefore used here. 

Details of each sample are given in Table 2, and cuts 
in color, magnitude and redshi ft are shown in Figure [TJ 
Following IWillmer et al.l ()2006D . we define red and blue 
galaxies using the observed color bimodality in DEEP2 
with the following tilted cut in color-magnitude space (in 
AB magnitudes): 

{U -B) = -0.032(AfB + 21.62) -I- 1.035. (1) 

We do not allow this color cut to evolve with redshift 
within the redshift range used here. The color samples 
are shown in Figure [U where the upper panels are color- 
magnitude diagrams for all DEEP2 galaxies with 0.9 < 
z < 1.0. The main red and blue samples are defined 
as those galaxies lying redward or blueward of the line 
defined in Eqn. 1, with Mb < —20 (dashed lines in the 
upper left panel of Figure [1]) and in the redshift range 
0.7 < z < 0.925 for red galaxies and 0.7 < z < 1.05 for 
blue galaxies. While the redshift ranges are not identical 
for the main red and blue samples, differences in their 
clustering properties are due to the color differences and 
not the redshift ranges, as we have checked that using the 
same redshift range for both does not change the results. 
Using similar redshift ranges also minimizes the effects of 
cosmic variance, as any density fluctuations will (to flrst 
order) be in common between the red and blue samples. 
We define green galaxies as those within ([/ — i?) = 0.1 
of the red/blue dividing line (shown as dot-dash lines in 
Figure [T|) in the redshift range 0.7 < z < 1.0. The green 
galaxy sample contains a subset of galaxies from each 
of the red and blue galaxy samples; it is not defined as 
a distinct set of objects, which would require removing 
green galaxies from the main red and blue samples. 

We further divide the red and blue samples into finer 
color bins: red galaxies are divided into redder and bluer 
halves at (U — B) — 1.21, while blue galaxies are di- 
vided into three color samples at {U ~ B) — 0.65 and 
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{U ~ B) — 0.79. These color cuts are shown in the up- 
per right panel of Figure [T] These finer color cuts are 
independent of magnitude and do not have the a tilt in 
the color-magnitude diagram, unlike the red/blue divid- 
ing line. Our results do not change significantly if we 
instead define all cuts to be dependent on both color and 
magnitude. The magnitude range we sample here is not 
wide, and the number of galaxies near these dividing lines 
is small enough that it does not significantly affect our 
results. 

Within the red and blue galaxy populations we also 
construct subsamples as a function of threshold luminos- 
ity. For the brighter luminosity samples we allow the up- 
per redshift limit to increase with increasing luminosity, 
so as to include more galaxies and reduce measurement 
errors. The lower panels of Figure [1] show Mb as a func- 
tion of redshift for red (left) and blue (right) galaxies. 
The dashed lines indicate the various luminosity sam- 
ples defined in Table 2. However, as illustrated in Figure 
2 of iGerke et all (|2007f l. the Rab = 24.1 DEEP2 tar- 
get selection limit defines a restframe color-dependent 
selection which changes as a function of redshift and 
which is not evident in Figure [TJ Figure 2 of lGerke et al.l 
(|2007f ) contains {U - B)-Mb colo r-magnitude diagrams 
of DEEP2 galaxies in narrow redshift bins of 5z = 0.05. 
The redshift-dependent selection is clearly evident, with 
the result being that red galaxies are not included in the 
survey at a given Mb which is complete for blue galax- 
ies. This results in a lower upper redshift limit of the red 
samples used here compared to the blue samples of the 
same luminosity, though the mean redshift of the samples 
does not vary by more than 5z = 0.2. 

Number densities for each sample are given in Table 
2; these are derive d from the DEEP2 luminosity func- 
tion split by color (jWillmer et al .11200 6*). To account for 
redshift failures we use 'minimal' weighting for the blue 
galaxies and 'average ' weighting for the red galaxies (see 
IWillmer et all (|2006h for details). We estimate the num- 
ber density for the green galaxy sample by multiplying 
the number densities for the main red and blue samples 
by the percentage of galaxies in each color that are in 
the green valley in the redshift range 0.75 < z < 0.90 
(12.2% of blue galaxies and 21.6% of red galaxies). The 
dominant error in the number density estimate derived 
from the luminosity funct ion is the erro r on (jf; using 
the errors quoted in Willme ir" et al.l (|2006f ), the fractional 
errors on the number densities for the main red and blue 
samples (with Mb < -20.0) are -1-23%/ - 5% for blue 
galaxies and -1-40%/ — 55% for red galaxies. 

2.3. Coadded Spectra 

Both locally and to at least z = 1 galaxies are found to 
be bimodal not on ly in terms of their r estframe color 
distrib ution (e.g.. IStrateva et al.l 120011 : iBlanton et al.l 
I2003at iBell et al.l [20041) but also in morphological and 
spectral type fe.g.. iMadgwick et all 2002, 2003b). Lo- 
cally, there are strong correlations between morphol- 
ogy, spectral type and restframe color, with some scat- 
ter that likely refiects both intrinsic variation and the 
difficulty of unequivocally determining morphological 
type. Defining galaxy samples by spectral type in 
DEE P2 is not as st raightforward as using restframe col- 
ors (Madgwick et al. 2003b ), due to the limited spectral 
range of the DEEP2 spectra, which is a consequence of 



working at high resolution. We therefore define galaxy 
samples here by restframe color and luminosity; how- 
ever, we show in Figure H] coadded rest-frame spectra 
for galaxies in each of our color samples in order to in- 
vestigate and illustrate the physical nature of the opti- 
cal emission for galaxies in each sample. These coadded 
spectra demonstrate that our color divisions correspond 
directly to classical spectral types. The coadditions cor- 
rect for the relative probability that each galaxy was se- 
lecteding for observation by DEEP2, but otherwise all 
objects are given equal weight. To maximize signal-to- 
noise and robustness, pixels affected by night sky emis- 
sion are deweighted and those affected by instrumental 
artifacts are removed. 

The top row of Figure |2| includes the coadded spec- 
tra for the main blue and red galaxy samples. The blue 
galaxies display strong emission lines, primarily from star 
formation ([Oil] A3727, H7, H/3, [OIII] AA4959, 5007), 
along with weak Ca H-l-K AA3934, 3968 absorption and 
stronger Balmer absorption features. The main red 
galaxy sample has a small am ount of [OH] an d [OIII] 
emission, primarily from AGN (lYan et al.ll2006l) . as well 
as Ca H-t-K absorption and weak Balmer absorption fea- 
tures. 

We also show coadded spectra for galaxy samples se- 
lected in finer color bins. Within the blue cloud, the 
bluer galaxies have much stronger star-forming emission 
lines in equivalent width, while the redder galaxies have 
lower EW lines; relative to the normalized continuum; in 
luminosity the redder of the blue galaxies have stronger 
[Oil]. 

The green galaxies, composed of both blue and red 
galaxies, have lower EW star-forming lines than the blue 
galaxies. The green galaxy spectrum does not appear 
to be a simple mix of the red and blue galaxy spec- 
tra, as the line ratios are different. The green galaxy 
spectrum has a much higher [0III]/II/3 ratio than the 
blue galaxy spectrum or the average of the 'blue:redder' 
and 'red:bluer' spectra; this indicates a higher fraction 
of AGN in the green population relative than in the blue 
cloud as a whole (Salim et al. 2007; Nandra et al. 2007), 
and also shows that the green galaxies are not a simple 
average of the redder of the blue cloud galaxies and the 
bluer of the red sequence galaxies. The high [OIII]/II/? 
ratio in the green valley population is likely dominated 
by the red galaxies in the sample (Weiner et al. 2007). 
Within the red sequence, the bluer galaxies have a higher 
[0III]/II/5 ratio, indicating more AGN activity, and the 
redder galaxies have stronger Ca H-l-K absorption fea- 
tures. 

3. METHODS 

3.1. Measuring the Two-Point Correlation Function 

The two-point correlation function ^(r) is defined as a 
measure of the excess probability dP (above that for an 
unclustered distribution) of finding a galaxy in a volume 
element dV at a separation r from another randomly- 
chosen galaxy, 

dP = n[l + ^ir)]dV, (2) 

where n is t he mean numb er density of the galaxy sample 
in question (jPeebles '' 19801. 

For each galaxy sample we construct a randomly- 
distributed catalog with the same overall sky cover- 
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age and redshift distribution as the data. We then 
measure the two-poin t correlation function using the 
iLandv fc Szalavl (fl99l estimator, 



1 



DD 



no 



2DR 



no 



RR 



(3) 



where DD, DR, and RR are counts of pairs of galax- 
ies (as a function of separation) in the data-data, data- 
random, and random-random catalogs, and no and hr 
are the mean number densities of galaxies in the data and 
random catalogs. This estimator has be en shown to per- 
form as well as the Hamilton estimator (jHamiltonI 1 1 993( 1 
but is preferred here as it is relatively insensitive to the 
size of the random catalog and handles edge corrections 
well ([Kcrscher, Szapudi, & Szalav 2000). 

To estimate the cross-correlation function between two 
galaxy samples, we measure the observed number of 
galaxies from a given sample around each galaxy in an- 
other sample as a function of distance, divided by the 
expected number of galaxies for a random distribution. 
For this we use the simple estimator 



DjD2 
D1R2 



1, 



(4) 



where D1D2 are pairs of galaxies between the two sam- 
ples and D1R2 are galaxy-random pairs between one 
galaxy sample and the random catalog of the other sam- 
ple, where the pair counts have been normalized by n^i 
and UR. 

The DEEP2 redshift success rate is >70% overah (de- 
fined as the percentage of galaxies targeted for spec- 
troscopy that have a well-determined redshift) and is not 
entirely uniform across the survey; some slitmasks are 
observed under better conditions than others and there- 
fore yield a slightly higher completeness. We only use 
regions of the survey with a redshift success rate >65%, 
and the spatially-varying success rate is taken into ac- 
count in the window function, which is applied to the 
random catalog to ensure that it has the same spatial 
distribution as the survey. The two-dimensional window 
function of the DEEP2 data in the plane of the sky re- 
flects the probability of observing a galaxy and takes into 
account the overall outline of the survey and the geom- 
etry of the overlapping slitmasks as well as vignetting 
in the DEIMOS camera and gaps between the DEIMOS 
CCDs. The varying redshift success is also taken into 
account such that regions of the sky with a higher com- 
pleteness have a correspondingly higher number of ran- 
dom points. This ensures that there is no bias introduced 
when computing correlation statistics. We also mask the 
regions of both the random and mock catalogs where the 
photometric data have saturated stars or CCD defects. 

Redshift-space distortions due to peculiar velocities 
along the line of sight significantly affect estimates of 
^(r). At small separations, random motions within a 
virialized overdensity cause an elongation along the line 
of sight ("fingers of god"), while on large scales, coherent 
infall of galaxies into potential wells causes an a pparent 
contr action of structure along the line-of-sight (|Kaiseij 
Il987t ) . While these distortions can be used to uncover in- 
formation about the underlying matter density and ther- 
mal motions of the galaxies, they complicate a measure- 
ment of the two-point correlation function in real space. 



In order to separate the effects of these redshift-space 
distortions and uncover the underlying real-space clus- 
tering properties, we measure ^ in two dimensions, both 
perpendicular to a nd along the line of sight. Following 
'Fishe r et all lfigg?). we define vectors vi and V2 to be the 
redshift-space positions of a pair of galaxies, s to be the 
redshift-space separation (vi — V2), and 1 = ^(vi -I- V2) 
to be the mean coordinate of the pair. We then define 
the separation between the two galaxies along (tt) and 
across (r^) the line of sight as 



si 

IT' 



(5) 



(6) 



To es timate ^(rp,7r), we apply the ILandv fc Szalavl 
(|1993f ) estimator to pair counts subdivided over a grid 



m rp and tt. 



3.2. Deriving the Real-Space Correlations 

While C(rp,7r) contains useful information about pe- 
culiar velocities, we also wish to measure the real-space 
correlation function, ^(r). To recover ^(r) we use a pro- 
jection of ^(^pjTt) along the rp axis. As redshift-space 
distortions affect only the line of sight component of 
£,{rp, tt), integrating over the tt direction leads to a statis- 
tic Wp{rp), whic h is independent of reds hift-space distor- 
tions. Following [Davis &: Peeblei (|1983( ). 

pQC poo 

Wp{rp) = 2 / dTT e(rp, tt) - 2 / dy ^[(r^ + y^/^ 
Jo Jq 

where y is the real-space separation along the line of 
sight. If^(r) is modelled as a power law, £_{r) = {r/ro)~^, 
then To and 7 can be readily extracted from the projected 
correlation function, Wp{rp), using an analytic solution to 
Equation [71 



Wp{rp) 



7r(i)r(i^) 



r(i) 



(8) 



where F is the usual gamma function. A power-law fit 
to Wp{rp) will then recover rg and 7 for the real-space 
correlation function, £,(r). 

In practice, ^(r) is not expected to be a power law at 
very large scales (>> 10 Mpc), nor can we measure 
(^(rp,7r) accurately to infinite separations as assumed in 
Equation[7l Here we integrate Wp{rp) to TTmax — 20 
Mpc, as £,{rp, n) becomes noisy at larger separations. It 
is not appropriate to then apply equation [8] directly; for 
power-law correlation functions, that equation becomes 
increasingly inaccurate at large rp (failing significantly 
where rp/TTmax ^ 0.25). Instead, we must compare the 
observed Wp{rp) to the integral of C(rp,7r) over separa- 
tions within TTniax, as predicted for a given rg and 7. 

This prediction is complicated by the presence of 
red shift-space dis torti ons due to coherent infall of galax- 
ies (iKaiseil (119871); seelHamiltonI (jl992l ) and section 4.1 of 
[Hawkins et al.l ([20031 ) for the relevant equations for cor- 
relation function analyses). If TTmax were infinite, these 
distortions would have no effect, as they merely change 
the line-of-sight separations of pairs of galaxies; but that 
is not the case here. The strength of these distortions 
depends on (3 = il^^/b, where ilm is defined at the mean 
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redshift of the sample (not z = 0) and b is the hnear bias 
between the clustering of galaxies and dark matter. 

We therefore recover ro and 7 as follows. We first 
measure Wp{rp) in the data by integrating ^(rp,7r) to 
TTmax = 20 Mpc. Wc then estimate Wp{rp) for dark 
matter particles at the mean redshift of the data sample , 
using the publicly-available code of ISmith et al.l ()2003D 
(not including redshift-space distortions), where we inte- 
grate the dark matter £,{rp,7T) to TTmax = 20 Mpc. 
We then estimate the linear galaxy bias from the ratio 
of these quantities: &^ = [wp]gai/[wp]daikmattor- We use 
the average value of this b over scales = 1 — 10 
Mpc and use ri,„ at the mean redshift of the sample (for 
an LCDM model with flm{0) = 0.3), to provide an esti- 
mate of the redshift-space distortion parameter /3. Using 
i^m(O) = 0.24 rather than 0.3 increases our recovered rg 
values by 1% and does not change 7. We can then pre- 
dict the observed Wp{rp) for any point in this grid by 
integrating C(rp,7r) to 7rmax=20 Mpc. The best-fit 
ro and 7 values are those that minimize the difference 
between the data and the model prediction. 

Although our estimated bias, and hence the /3 used in 
this procedure, may not perfectly match the bias corre- 
sponding to the best-fit values of ro and 7, the measured 
values prove to be insensitive to that estimate. Increas- 
ing the bias by 20% increases rp by 1% and decreases 7 
by 1%. This method for recovering ro and 7 assumes that 
^(r) is a power law only to a scale of TTmax and results in 
ro and 7 values within a few percent of those obtained 
using Equation [U for our value of TTmax = 20 Mpc. 
Deviations from Equation [8] are significant only on larger 
scales, where rp/n^ax ^ 0.25. We also consider an alter- 
native method for recovering correlation functions in the 
next section. 

Errors on Wp{rp) are calculated using the standard er- 
ror across the 10 separate data pointings. We do not 
use a full covariance matrix in the fits. Errors on ro and 
7 are derived from jacknife resampling of the separate 
pointings and therefore take into account the covariance 
among the Vp bins. 

3.3. Comparison with an Alternative Clustering 
Statistic 

Given that it becomes necessary to model the ef- 
fects of the truncation of Wp and include redshift-space 
distortions when comparing theory with observations, 
Padmanabhan, White and Eisenstein (2007, hereafter 
PWE07) propose an alternative statistic, uj{Rs), which 
has better convergence properties to the real-space quan- 
tities. This statistic is defined as 

uj{Rs) = 2Tr J rpdrp G{rp, Rs)wp{rp) 

= 27r / rpdvp G{rp,Rs) I d7r^(rp,7r) (9) 

where again Vp and tt are transverse and line of sight 
directions, respectively. The filter function, G{rp,Rs) is 
chosen to be non-zero only for Rg < 1 and to have zero 
integral (when integrated against an r^ measure). We 
follow PWE07 and adopt G(rp) = Rj^G{x = rp/R^) 
with 

G{x) = a;4(l - a;2)2 ( 1 - x"^) x < 1 (10) 
= x > 1 . (11) 



This statistic has the further advantages of being un- 
binned and having a well- localized kernel in ^(r), prob- 
ing the correlation function at a scale ~ Rg. We note 
that w(i?s) can be thought of as performing the Abel in- 
version from Wp to ^, while avoiding differentiating noisy 
data. 

In order to test the robustness of the methods pre- 
sented above, below in Section 4.1 we compute w for the 
main red and blue galaxy samples to compare with the 
results obtained using the methods in S ection 3.2. To 
compute uj we use the iLandv fc Szalavl (il993.) estima- 
tor for the two-dimensional correlation function, and the 
techniques in PWE07 to rewrite the integrals in Eqs.[9]as 
weighted pair sums; the line of sight integral is truncated 
at TTmax = 40 Mpc. Uncertainties in lu are estimated 
by jackknife resampling the ten DEEP2 pointings. Since 
UJ is unbinned, one can choose arbitrarily finely spaced 
bins, with the disadvantage that adjacent bins become 
highly correlated. We have verified explicitly that our 
fits for ro and 7 given below are insensitive to the par- 
ticular choice of binning we adopt. 

It is straightforward to relate ui to the three- 
dimensional correlation function, ^. For a power-law 
correlation function, ^(r) = (r/ro)~^, and our choice of 
G{rp, Rs), this relation becomes 

^ 2^3/2r(7_l) 4(^_^) /^x - 

r(l) (7-7)(9-7)(ll-7)(13-7) Uoy 

(12) 

3.4. Systematic Biases due to Slitmask Design 

When observing with multi-object slitmasks, the spec- 
tra of targets cannot be allowed to overlap on the CCD 
array; therefore, objects that lie near each other in the 
direction on the sky that is perpendicular to the wave- 
length direction on the CCD cannot be simultaneously 
observed. This results in undersampling the regions of 
the sky with the highest density of targets. To reduce 
the impact of this bias, adjacent DEEP2 slitmasks are 
positioned approximately a half-mask width apart, giv- 
ing each galaxy two chances to appear on a mask; we also 
adaptively tile the slitmask locations to hold constant 
the number of targets per mask. In spite of these steps, 
the probability that a target is selected for spectroscopy 
is diminished by ~ 25% if the distance to its second 
nearest neighbor is less than 10 arc seconds (for techn i- 
cal details of the slitmask d esign see lDavis et al.l ()2003D : 
iDavis. Gerke. fc NewmanI (|2005l )). This introduces a 
roughly predictable systematic bias which leads to un- 
derestimating the correlation strength on small scales. 
The objects that are in conflict are often not at the same 
redshift, however, such that the effect is not particularly 
large. 

To correct for thi s effect we use the mock catalogs of 
lYan. White, fc Coil' (20041) with galaxy colors added as 
described in [Gerkc ct aLr (l2007l) . We measure the ra- 
tio of both ^(rp, tt) and Wp{rp) in the mock catalogs be- 
tween samples of galaxies from catalogs with and without 
the slitmask target selection algorithm applied. We then 
multiply the measured ^(r^, tt) and Wp{rp) in the data by 
this ratio, which is a smooth function of scale and varies 
from ~ 20% on the smallest scales to 2% on scales greater 
than rp = 1 Mpc. We tested mock catalog samples 
with identical redshift, luminosity and color ranges as 
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each data sample used here and found no significant dif- 
ferences in the multipHcative corrections as a function 
of color or luminosity, such that the same correction de- 
rived from the largest sample is our best estimate for all 
other samples. The effect of this slitmask correction is 
small; it increases tq by 1.5% and 7 by 2.5%. We include 
an additional error due to this correction, which is added 
in quadrature to the error on Wp{rp). This error is esti- 
mated from the variance on the correction in the mock 
catalogs and is a function of scale, varying from 10% at 
rp = 0.05 h"^ Mpc to 2% on scales Vp > 1 h^^ Mpc. 

4. GALAXY CLUSTERING RESULTS 
4.1. Blue and Red Galaxies 

We first compare the main blue and red galaxy sam- 
ples, both of which have Mb < —20. The top left panels 
in Figure [3] show ^(rp,Tr) for these two samples, with 
contours indicating constant probability, where the dark 
line is ^ = 1. There are several clear trends present in 
these diagrams; red galaxies are seen to be more clus- 
tered than blue galaxies (a given ^ is seen at larger sep- 
arations) and show much stronger fingers of god, as seen 
in the elongation along the tt direction at small Vp sepa- 
rations (<1 Mpc). Both of these results reflect the 
fact that a color-density relation is in place at ^ ~ 1, with 
red galaxies residing more in overdense regions such as 
groups and clusters than blue galaxies. Infall on large 
scales is plainly seen in these ^(r^, tt) diagrams for blue 
galaxies, as indicated by the flattening of the contours 
along the line of sight. Coherent infall may be present 
for red galaxies as well, though it is harder to detect due 
to the larger fingers of god (see Section 4.9 for more de- 
tails). Redshift-space distortions are significant out to 
TT = 20 Mpc, such that it is important to not use 
the usual power-law estimator (Equation [8]) to infer tq 
and 7 from uupirp), as this equation assumes that one has 
integrated ^(rp, tt) to tt = 00. 

The top left panel of Figure [3] shows the projected cor- 
relation function, Wp{rp), for the main red and blue sam- 
ples. Red galaxies exhibit a significantly steeper slope 
than blue galaxies, and both samples are well approx- 
imated by a power law within the error bars. Power- 
law fits are given in Table 2, fitting the data on scales 
rp = 0.1 - 20 h-^ Mpc; the smallest Vp bin (0.05 h'^ 
Mpc) is not used in these fits due to the rise on small 
scales which is seen in some samples (see Section 4.3 for 
further discussion of deviations on small scales). Includ- 
ing the smallest bin in these fits, or fitting only on scales 
larger than rp = 1 Mpc, does not change the results 
for the main blue and red galaxy samples. 

The correlation length for red galaxies is found to be 
significantly higher than for blue galaxies: tq = 5.25 ± 
0.26 h-^ Mpc compared to ro = 3.87 ± 0.12 h"^ Mpc. 
These values for rp at z ~ 1 are somewhat lower than 
what is found for local samples (Table 1); however, tq has 
only increased by ~5-15% by z ^ 0.1 for blue galaxies 
and ~ 10-20% for red galaxies. Furthermore, there has 
been very little evolution in the slope of the correlation 
function for either red or blue galaxies since z ^ 1. We 
discuss the implications of this in Section 6. 

Our ro values for blue and re d galaxies a re somewhat 
higher than those r eported by ICoil et al] ((2004b ) and 
iMeneux et al.l ()2006[ ). though generally within the large 



errors in those findings. ICoil et al] ()2004bD use early data 
from a portion of one DEEP2 field and iMeneux et al.l 
(I2OO6I) use the VVDS Deep Survey first epoch data; both 
of these sa mples are roughly 1/ 6 of the dataset used here, 
and neither ICoil et all (pOOib) nor IMeneux et all (|2006l ) 
constructs volume-limited subsamples. For blue galax- 
ies, our rn value found here is 2a and 5cr higher than 
ICoil et al.l (|2004bD and IMeneux et all (|2006D . while for 
red galaxies our ro value is la and 2a higher, respec- 
tively, where the er rors have been added in quadrature. 
iBrown et al l (120031) measure the clustering of red galax- 
ies in the NO AO Deep Wide-Field Survey, using photo- 
metric redshifts. They find at z = 0.85 that ro = 6.7±0.8 
for red galaxies with —21.5 < Mn < —2 0.5 (roughly 
L/L->f) , which is 2a higher than our result. iHeinis et al.l 
measure the clustering of UV-selected galaxies, 
again using photometric redshifts, and find that at 2: ~ 
0.9 ro = 4.92 ± 0.5 h'^ Mpc and 7 = 1.7 ± 0.09. Their 
clustering scale-length is 2a higher than what we find 
here for star-forming galaxies. 

We also use the estimator of PWE07, as described in 
Section 3.3, to measure w for the main blue and red 
samples; the results are shown in Fig. |5l We have not 
attempted to correct for slitmask effects in the u) esti- 
mator; these effects are small, as discussed above. Using 
Eqn. 1121 we fit for ro and 7 for both galaxy samples. The 
resulting power-law fits are plotted in Fig. |5| (solid lines), 
and are compared with the fits obtained from fitting to 
Wp (dashed lines). The two fits are clearly consistent 
with each other, given the measurement errors. For the 
blue galaxy sample, we obtain ro — 4.45 ± 0.54 Mpc 
and 7 — 1.47 ± 0.11, while for the red galaxies, we find 
ro = 5.94 ± 0.57 h'^ Mpc and 7 = 1.93 ± 0.08. These 
ro and 7 values are consistent with the estimates derived 
above from fitting to Wp. In fact, because the correla- 
tion coefficient between ro and 7 is large and negative 
—0.8), difference between the results of the two esti- 
mators are much less significant than the projected errors 
imply. The error bars on ro and 7 are larger using the 
uj estimator, which is not surprising given that lo is an 
inversion of Wp to a localized integral of ^(r). 

4.2. Red-Blue Cross-Correlation 

In addition to measuring the auto-correlation of blue 
and red galaxies separately, the cross-correlation between 
the two samples can be used to further understand the 
spatial relationship between blue and red galaxies. The 
cross-correlation between two samples measures the clus- 
tering of one type of object around the other and pro- 
vides information on the mixing of the populations. On 
large, linear scales, in the two-halo regime, blue and red 
galaxies should be well-mixed and trace the same over- 
all large-scale structure, such that the cross-correlation 
should equal the geometric mean of the auto-correlations 
of each sample. On small scales, however, within the 
one-halo term, deviations from the geometric mean en- 
code information about the differences in the halos the 
two samples populate. For example, if some halos con- 
tain almost exclusively red galaxies and others only blue 
galaxies, then the cross-correlation function would fall 
below the geometric mean. The blue-red galaxy cross- 
correlation function can therefore be used to distinguish 
between different star formation quenching scenarios, as 
it traces the distribution of quiescent and actively star- 
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forming galaxies around each other. 

Here we measure the cross-correlation of the main 
blue and red galaxy samples to investigate in particu- 
lar the mixing of blue and red galaxies on small, one- 
halo-dominated, scales {vp < 1 Mpc). The two- 
dimensional ^(rp,7r) cross-correlation function is shown 
in the lower left panel of Figure [3] and the projected cor- 
relation function, Wp{rp), is shown in the upper right 
panel of Figure 21 The cross-correlation is intermediate 
between the auto-correlation of blue and red galaxies. 
The best fit scale length is rg = 4.37 ± 0.41 Mpc 
and slope is 7 = 1.83 ± 0.04. Fingers of god are clearly 
seen in C(rp, tt), though the coherent infall on large scales 
appears to be more similar to what is seen for the red 
sample than the blue sample. 

Figure [6] shows Wp(rp) for the main blue and red galaxy 
samples and the red-blue cross-correlation function di- 
vided by a reference power law Wp{rp) with tq = 4.0 
Mpc and 7 = 1.8. Also plotted are deviations from this 
power law for the geometric mean: 

Wp(Tp). red— blue — \j Wp{rp)j,edWp {rp)blue- (13) 

While the cross-correlation (dot-dash line) is consistent 
with the geometric mean (dotted line) within the error 
bars, on scales rp < 1 Mpc (the one-halo term) the 
cross-correlation is closer to the blue Wp{rp) than the red. 
This reflects a deficit of blue galaxies near red galax- 
ies, which could be due to the centers of groups being 
occupied preferentially by red galaxies. This is consis- 
tent with the observe d cross-correlation between DEEP2 
groups and galaxies (jCoil et al.l l2006b\ which is higher 
for red galaxies than blue galaxies, again indicating that 
red galaxies tend to be at the centers of groups. 

A similar effect has been seen at z = 0.1 in the SPSS 
red-bl ue cr oss-correlation function. Both IZehavi et al.l 
(12003) and ISwanson et~an (|2007f ) find that the red-blue 
cross-correlation amplitude is lower than the geometric 
mean of th e red and blue clustering on small scales. 
iWang et all ((200I measure the cross-correlation as a 
function of luminosity and find a suppression relative to 
the geometric mean for faint red and fai nt blue g alaxies , 
but not for brighter samples. Weinma nn et al.l (|2006bf ) 
find a similar effect in their SDSS group catalog, which 
they term 'galactic conformity', in which the colors of 
satellite galaxies are correlated with the colors of their 
central galaxies; i.e., red satellite galaxies tend to be 
around red central galaxies, and the same is seen for blue 
galaxies. 

4.3. Luminosity Dependence for Blue and Red Galaxies 

The luminosity dependence of clustering i n the over- 
all DE EP 2 galaxy sample was presented in ICoil et all 
(|2006ar ): here we expand upon that study by comput- 
ing the luminosity dependence of clustering within the 
blue and red galaxy populations separately. The median 
Mb of each sample is given in Table 2. For compari- 
son, &t z = 0.9, M* is Mb = -20.54(-H0.03/ - 0.01) for 
blue galaxies and Mb = — 20.35(±0.03) for red galax- 
ies (from Willmer et al. 2006, where we use the 'mini- 
mal' weighting for blue galaxies and 'average' weighting 
for red galaxies). The main blue sample has a median 
Mb value near M* , while the blue luminosity subsam- 
ples range from 0.44 magnitudes fainter than M* to 0.74 



magnitudes brighter than M* . The main red sample is 
0.35 magnitudes brighter than M* (in the median) for 
red galaxies, and the red luminosity subsamples range 
from ~ M* to 1 magnitude brighter than M* . 

Figures iTl and IHl show ^(rp,7r) for four luminosity 
threshold samples for blue and red galaxies. Within the 
blue population, there is almost no difference between 
the faintest two samples, while the brightest samples may 
show somewhat stronger fingers of god. All of the blue 
luminosity samples show coherent infall on large scales. 
The red ^(rp, tt) shows almost no difference between any 
of the bins; the brightest is noisy, however, due to the 
small sample size. Any correlation between the strength 
of the small-scale redshift-space distortion and luminos- 
ity within the red galaxy population, for the luminosity 
range that we p robe here, is too weak for us to detect. 
iLi et all (|2006b[ ) find that the pairwise velocity disper- 
sion (see Section 4.8), which quantifies the amplitude 
of the fingers of god, increases with luminosity for blue 
galaxies in SDSS, while red galaxies do not show any lu- 
minosity dependence, qualitatively similar to our results 
at z ^ 1. 

The projected correlation function, Wp{rp), is shown 
in Figure 1^1 for blue subsamples (top row) and red sub- 
samples (bottom row) as a function of luminosity. Again 
there is little difference in clustering between the fainter 
blue samples (upper left panel); however Wp{rp) displays 
a significant rise on small scales for the brighter blue sam- 
ples (upper right panel). The red samples show no cor- 
responding rise on small scales for the brighter samples 
(lower right panel). Deviations from a power law for the 
brighter samples are shown in Figure ITOl on scales Vp < 
0.2 Mpc there is a clear rise in the correlation func- 
tion of blue galaxies. A similar incr ease on sm a ll scale s 
was seen for ah DEEP2 galaxies in ICoil et all (|2006aD . 
This is likely due to brighter galaxies having, on average, 
a higher probabilit y of being a central galaxy rather than 
a sat e Uite galaxy (IZehavi et al."2 005tlMandelbaum et al. 
20061 Ivan den Bosch et al. 2007: IZheng. Coil, fc Zehavi 



20071 ) , which leads to a relative rise in the one- halo term 



compared to the two-halo term (Tinker et al, in prep.). 
There may also be a contribution from interaction -driven 
starbursts fe.g.. iBarton. Geller. Sz KenvonI I2OOOI ) . The 
lack of a strong rise on small scales in the brighter red 
samples likely results from red DEEP2 galaxies having 
a higher satellite fraction than blue galaxies, at a given 
luminosity. We note that while the brighter blue galax- 
ies are more likely to be central galaxies, they are not 
necessarily in more massive halos, as shown by the lack 
of significant dependence of tq on luminosity within the 
blue cloud. 

Power-law fits are given in Table 2 where, as above, the 
scales used to fit a power law are Vp = 0.1 — 20 h"-^ Mpc. 
Fitting on scales rp = 1 — 20 Mpc instead does not 
significantly change the results. The results are shown 
in the upper panels of Figure 1111 There is no significant 
dependence in the correlation length, tq, or slope, 7, on 
luminosity for either blue or red galaxies, within the error 
bars. 

ICooper et al.l (|2006f ) use local galaxy overdensity mea- 
sures to determine the environment of DEEP 2 galaxies 
and conclude that within the luminosity range sampled 
by DEEP2 there is no significant dependence on luminos- 
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ity for red galaxies, although there is for blue galaxies, in 
that brighter blue galaxies are in more overdense regions. 
Within the luminosity range pr obed here, the relat ive 
bias of blue galaxies is found by ICooper et al.l ()2006D to 
increase by ~30%, which is steeper than the trend found 
here using the correlation function. The brighter blue 
samples do show an increase in the small-scale cluster- 
ing amplitude relative to the fainter blue samples, how- 
ever, which is not re flected in the measured rg values. 
ICucciati et al.l (2006) also measure local overdensity at 
z ~ 0.5 — 1.5 in VVDS data and find that the fraction of 
red galaxies at a given overdensity depends significantly 
on luminosity. We do not find such a trend here with our 
larger dataset. 

The lower panels of Figure [TT] compare the clustering 
scale-length of red and blue DEEP 2 gala xies as a func- 
tion of LjL* with local results from 2dF (|Norberg et al.l 
l2002bH for passive and active galaxies, defined usin g spec- 
tral types, and from SDSS (jZehavi et al.l 120051 ) . using 
restframe colors. For this comparison we use = 
-20.67 for aU galaxies at z = 0.9 ( Willmer et al.ll2006f) 
and Ml, = -19.66 and Mr = -20.44 for all galax- 

ies at z = 0.1 (iNorberg et al.l l2002at 'Blanto n et~al 
l2003bD . At z = 1, IZfi havietalJ feoOS); Li et al. 
jlOOl^); IWang et all (|2007f l all detect significant lumi- 
nosity depende nce in the clustering o f both blue and red 
galaxies, while ISwanson et al.l (l2007l) detect luminosity 
dependence for red galaxies onlv. INorberg et"an (j2002a l 
detect significant luminosity dependence in the clustering 
of both passive and active galaxies as well. 

This luminosity dependence for local blue and red 
galaxies is significant at i > L* but not at L ^ L*; 
there is also a significant increase in the clustering for 
faint red galaxies, but that is below the luminosity range 
that we probe here. The DEEP2 sample docs not probe 
L > 2L*, such that we can not test the luminosity depen- 
dence at large L/L*. Within the range that we do probe, 
it does not appear that either blue or red galaxies have 
a detectable luminosity dependence in their autocorrela- 
tion properties. This is consistent with local results in 
the same L/L* range. 

Our result that the clustering amplitude does not 
strongly dep end on lumino s ity fo r red galaxies at 2 ~ 1 is 
at odds with iBrown et all (j2003f ). who find a significant 
luminosity dependence in the clustering of red galaxies 
at z = 0.3 — 0.9, for samples between Mr = —20 and 
Mr = —22 (in Vega mag nitudes), using pho tometric 
redshift data. for the iBrown et all (|2003f) sample 

is ~ —21.0 (M. Brown, private communication). Their 
value of ro at z ^ 0.6 for L/L* ~ 1 red galaxies is 
6.3 ± 0.5 Mpc, which is only somewhat higher than 
that foun d here (tq = 5.25 ± 0.26 at z = 0.82) or at 
z - 0.1 bv lZehavi etaJ] (|2005h (ro = 5.67 ±0.37). How- 
ever, they find a steep trend in tq with M^, with tq 
ranging from 4.4 ± 0.4 h'^ Mpc to 11.2 ± 1.0 h'^ Mpc 
for L/L* ~ 0.4 — 2.5, which is not found here or lo- 
cally. No value of would lead to a consistent trend 
of rp with luminos ity when compared wit h our results or 
local re sults from IZehavi et al.l ()2005l ) or INorberg et al.l 
()2002bl ) . Estimati ng rp from an angular correlation func- 
tion, as done in Brown et al.l (f2003), requires very ac- 
curate knowledge of the redshift distribution; this may 
account for the discrepancy with the results from spec- 



tr oscopic surveys. 

ICoil et all (|2006af ) find that within the fuU DEEP2 
galaxy sample the clustering strength depends strongly 
on luminosity. This can still occur even if there is no lu- 
minosity dependence for red or blue samples separately, 
as the red galaxy fraction is a function of luminosity; the 
higher prevalence of red galaxies in the brighter samples 
leads to ro increasing as AIb decreases. 

4.4. Green Galaxies 

The clear bimodality both locally and at z = 1 in 
galaxy properties, whether defined by restframe color, 
morphology or spectral type, raises the question of how 
these two main types arise and what, if any, evolutionary 
connection exists between them. The observed buildup of 
galaxies on the red sequence since z = I (e . g . , iBell et aP 
12004 iFaber et al.ll2007l : IBrown et al.ll2007D suggests that 
blue galaxies are moving to the red sequence with time, 
though the details of this transition are not clear. The 
general outline, however, is that star formation begins to 
shut down or be quenched in blue galaxies, after which 
they passively evolve onto the red sequence. The star for- 
mation quenching process (es) involved may be causally 
related to the environment of th e galaxy (e.g. , ram 
pressure or tidal stripp ing of gas (jGunn fc CxottJ [l97l 
iBvrd fc Valtoneiil Il99(l )). such that galaxies in higher 
density regions end their star forming phase more read- 
ily, thus creating a color-density relation. Alternatively, 
the process may be inherent to the galaxy and depend 
on its age or stellar or halo mass, in which case the pro- 
cess may be correlated with environment but not caused 
by it (|Bundv et al.ll2006HCooper et al.ll2007a^ . To eluci- 
date what is causing this change, we can investigate the 
clustering properties of galaxies located at the transition 
region between these populations. Here we study 'green' 
galaxies located near the minimum of the observed color 
bimodality as defined in §2.2. 

The green galaxy ^(rp,7r), shown in the upper right 
panel of Figure|3l exhibits an intermediate clustering am- 
plitude between the blue and red galaxy samples, as well 
as intermediate-strength fingers of god. The green pop- 
ulation displays somewhat similar infall on large scales 
as blue galaxies, though the contours are noisy on large 
scales due to the relatively small sample size. 

The projected correlation function for green galaxies is 
shown in the lower left panel of Figure (H plotted with 
Wp{rp) for red and blue galaxies for comparison. On 
large scales {vp > 1 h^^ Mpc) the green galaxies have 
a clustering amplitude similar to that of red galaxies, 
while on small scales the amplitude is closer to that of 
blue galaxies. The large-scale agreement with the red 
galaxy population implies that green galaxies reside in 
or on the outskirts of halos of similar mass as red galax- 
ies, as they have the same clustering for the two-halo 
term. Green galaxies are therefore generally in the same 
overdense regions as red galaxies. However, the one-halo 
term for green galaxies is lower than for red galaxies, 
which likely reflects a lower radial concentration of green 
galaxies within their parent halos. The overall picture is 
consistent with galaxies having green colors while on the 
outskirts of halos and becoming redder as they reach the 
center of the overdensity. These results point towards 
green galaxies being a distinct population and not a sim- 
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pie mix of red and blue galaxies. 

To investigate these trends further, we compute the 
cross-correlation between the green galaxy sample and 
the main blue and red galaxy samples. The ^(?'p,7r) 
diagrams (lower right panels of Figure [3|) show more 
significant infall on large scales in the green-blue cross- 
correlation than the green-red cross-correlation (compar- 
ing, for example, the ^ = 1 contours). The green-red 
cross-correlation is very similar to the blue-red cross- 
correlation (though with smaller fingers of god), indi- 
cating that green galaxies have comparable kinemat- 
ics to blue galaxies. There are also smaller fingers of 
god in the green-red cross-correlation than in the red 
auto-correlation. These results are consistent with green 
galaxies residing in the same overdense regions as red 
galaxies, on average, but with a lower probability of 
being at the cores of those over densities, compared to 
red galaxies. This would result in green galaxies hav- 
ing smaller fingers of god than red galaxies, and a lower 
auto-correlation on small scales. It is possible that green 
galaxies reside on the outskirts of the same halos that 
red galaxies occupy but are still falling in towards the 
centers of the halos. 

The bottom right panel of Figure[3]shows the projected 
cross-correlation functions between green and blue or red 
galaxies. The green-red galaxy clustering amplitude is 
higher than the green-blue galaxy amplitude on all scales. 
On scales rp = 1 — 10 Mpc the green-red sample is 
39 ± 5% more clustered. On small scales, dominated by 
the one-halo term, the green-red amplitude is a few times 
higher than the green-blue amplitude: a factor of 2 ± 0.2 
times higher on scales ~ 0.5 Mpc and a factor of 
2.8 ± 0.4 times higher on scales ~ 0.25 Mpc. One 
would therefore be much more likely to find a red close 
neighbor (or satellite) galaxy to a green galaxy than a 
blue neighbor, and vice versa, if red and blue galaxy 
samples were of equal size. Another way to frame these 
results is that a green galaxy is less likely to be a central 
galaxy than a red galaxy of similar luminosity, but more 
likely than a comparable blue galaxy. 

We note that the green galaxy Wp{rp) likely can not 
be fit with an HOD model that assumes that the radial 
distribution of galaxies follows that of dark matter parti- 
cles or subhalos. The green galaxy radial distribution is 
clearly different from that of either red or blue galaxies, 
and most HOD models do not allow the radial distribu- 
tion to vary for different galaxy types. 

4.5. Finer Color Bins 

We next divide the blue and red galaxy samples into 
finer color bins to investigate the strength of clustering 
as a function of color within each population. Figure 
[T^ shows To and 7 as a function of color for both the 
main blue and red samples (triangles), the green sample, 
and for finer color bins. As explained in Section 2.2, the 
finer color bins are not all independent; this is shown in 
the figure with a dotted horizontal line showing the color 
range for each point (for clarity the color ranges for the 
main blue and red samples are omitted). In particular, 
the green galaxy sample contains subsamples of both blue 
and red galaxies. 

Within the red population (where we have two inde- 
pendent samples) neither tq or 7 depend significantly on 
color; across the red sequence we detect no change in 



the clustering properties of galaxies. Within the blue 
population, however, there is a strong dependence; red- 
der galaxies are more clustered than the bluer galaxies. 
There is no corresponding trend seen in the slope of the 
correlation function. The abrupt change in slope between 
blue and red galaxies likely reflects the different radial 
concentrations and satellite fractions of the two popula- 
tions and may indicate that star formation quenching is 
an effect associated with the halo properties of galaxies. 

Interestingly, the green galaxy tq is found to be the 
same as the red galaxy tq (within 1 a), while the green 
galaxy 7 is the same as for blue galaxies, not red galaxies. 
This is consistent with the interpretation in the previous 
section that green galaxies reside in similar mass halos 
as red galaxies, and are found generally in overdense re- 
gions, but are less concentrated within their host halos 
than red galaxies. They apparently do not populate the 
centers of overdensities. 

As discussed earlier, luminosity function studies find 
that a significant fraction of the progenitors of local 
galaxies on the red sequence must have been in the 
blue cloud a± z ^ 1 (|Bell et al.ll2"00l iBrown et al']|2007t 
iFaber et al.ll2007l ). Our clustering results imply that the 
galaxies that migrate to the red sequence from z = 1 
to z = are amongst the redder of the blue galaxies, 
and the green galaxies at z = 1, which are already as 
clustered as galaxies on the red sequence. 

These trends of finding a color-density relation within 
the blue cloud but not within the red sequence are similar 
to what has been locally using galaxy environment mea- 
sures (e.g., Hogg et al. 2003, 2004). There are no other 
published z ~ 1 clustering results with such fine color 
bins to compare to, and we defe r a compariso n with the 
DEEP2 environment results of iCooper et al.l ()2006D to 
Section 4.7. 

4.6. Galaxy Bias and Relative Bias Between Samples 

To facilitate comparisons to published clustering re- 
sults measured at other redshifts or with different statis- 
tics than the two-point correlation function, we calcu- 
late the relative bias between various galaxy color sam- 
ples. We define the relative bias as the square root of 
the ratio of Wp{rp) for two samples. As the relative bias 
is scale-dependent, we calculate the mean relative bias 
over two scales, rp = 0.1 — 15 h~^ Mpc ('all scales') and 
— 1 — 15 h~^ Mpc ('large scales'), which have mean 
scales of rp — 4.4 h^^ Mpc and rp — 6.5 h^^ Mpc, re- 
spectively. 

The relative bias of red to blue galaxies is brei — 
1.28 ± 0.09 on large scales and Kei = 1-44 ± 0.07 on 
all scales, reflecting the rise in the relative bias on small 
scales, Tp < 1 Mpc. We do not find any depen- 
dence on luminosity for the relative bias between red 
and blue galaxies. T his is in contrast to the results of 
iCucciati et al.l ()2006[ ) , who measure the color-density re- 
lation at z ~ 0.5 — 1.5 in VVDS data and find that the 
color-density relation steepens for brighter galaxies. No 
such trend is detected here. 

The lack of luminosity dependence in the relative bias 
found here could indicate that the ratio of the satellite 
fraction in red versu s blue galaxy s ample s does not de- 
pend on luminosity. iBerlind et al.l (|2005!) show that for 
central galaxies there is no color-density relation; this 
relation only exists for satellite galaxies. This implies 
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that the relative bias of red to blue galaxies should be 
a function of luminosity if the central to satellite galaxy 
fraction changes with luminosity differentially for red and 
blue galaxies. The fact that the data show no dependence 
in the relative bias on luminosity could imply that the 
ratio of the satellite fraction in red versus blue galaxy 
samples is fixed and does not depend on luminosity. 

The relative bias of red to green galaxies is brei = 
1.00 ± 010 on large scales and b^ei — 1-12 ± 0.11 on all 
scales. Given that red and green galaxies have different 
correlation slopes and redshift-space distortions, the rel- 
ative bias does not fully reflect the differences between 
these populations. It does, however, indicate that when 
two-halo-dominated scales are included, on average the 
clustering amplitude of green galaxies is perfectly consis- 
tent with that of red galaxies and is not lower. The 
relative bias of green to blue galaxies is found to be 
brei = 1-30 ± 0.11 on both large and all scales, with no 
scale-dependence. 

The relative bias found here between red and blue 
galaxies at z ~ 1 is very similar to what is mea- 
sured at lower redsh ift (see Table 1; brei is/ ^ — 10 
% lower at z - 1). iMadgwick et all (|2003af ) compare 
the clustering of early and late type galaxies (defined 
by spectral type) and find a relative bias of brei = 
1.45 lb 0.14 integrated to scales of r = 8 Mpc, while 
IWillmer. da Costa, fc Pellegrinil (|1998f ) fi nd a relative 
bias of ~ 1.4 between red and blue galaxies. IZehavi et al.l 
(|2002) find using SDSS data that the relative bias of red 
to blue galaxies is ~ 1.6, but is highly scale-dependent. 
This agreement with lower redshifts implies that the de- 
pendence of color on density at z = 1 is as steep as it 
is locally; the physical mechanisms responsible for the 
color-density relation are just as effective before z = 1 as 
they are after. The color-density relation was therefore 
established before z = 1 and subsequent changes may be 
absolute (in terms of the zeropoint of the relation) but 
not relative (which wou ld change the s lope) , within the 
errors. As discussed in ICooper et al.l (|2006D . given the 
low abundance of rich clusters at z = 1, the color-density 
relation seen in the DEEP2 data can not be the result 
of cluster-specific physics such as ra m pressure str i pping 
of gas or galaxy haras sment (e.g., iGunn fc GotH \197% 
IBvrd fc ValtoiienlfTQQOh . 

For a given cosmology we can also compute the absolute 
galaxy bias relative to the underlying dark matter den- 
sity field using results fr om N-body s i mulat ions. We use 
the power spectrum of iSmith et al.l (|2003D to estimate 
the dark matter clustering amplitude at the same red- 
shift as the data, for a cosmology with flm = 0.3, JIa = 
0.7, CTg = 0.9, and F ~ 0.21. This code is an analytic 
fit motivated by the halo model; we include a 5% error 
on the dark matter Wp(rp) at each scale in our calcula- 
tion of the absolute bias to reflect the uncertainty in the 
fit. The resulting absolute bias values averaged on scales 
= 1 — 10 Mpc are listed for each color and magni- 
tude sample in Table 2. For as = 0.8 the large-scale bias 
is 13% higher, while for erg = 1-0 it is 10% lower. Note 
that the absolute bias and relative bias are measured on 
different scales. 

We find that the absolute bias is b = 1.65 ± 0.15 
for red galaxies and b = 1.28 ± 0.04 for blue galax- 
ies, both with Mb < —20.0. This is consistent with 
iMarinoni et al.l (|2005[ ). who find that red galaxies with 



Mb < —20.0 at z ~ 0.8 have an absolute bias of 
b = 1.5 ± 0.6 and blue galaxies have b ~ 1.1 ± 0.6, 
measured on scales of r = 8 Mpc. The abso- 
lute bias of galaxies at z ^ 1 measured here is higher 
than at z = (e.g., IVerde et aT 2002; ' Hoekstra et al.l 
120021: iSeliak et al.ll2005l : iSimon et al. 2007i) . where galax- 
ies near L* are found to have a bias of 6 < 1.0. 

4.7. Comparison to Environment Studies 

To compare our results on the color dependence of 
g a.laxy clustering with the DEEP2 environment results 
of lCooper etldl (f200ll . we show in Figure [T3] the relative 
bias of each color sample used here, relative to the main 
blue sample, measured on scales Tp = 1 — 5 Mpc. 
The main red and blue samples are shown as triangles, 
while the finer color bins are shown as diamonds. Crosses 
with dotted error bars are the relative galaxy overdensity 
as measured by the en vironment statistic St, (described 
in iCooper et al.l ()2006D ) for galaxies with z = 0.75 - 1.0 
and Mb < —20.0, as a function of color, again normal- 
ized to the main blue sample. Note that this redshift and 
m agnitude rang e are n ot exactly the same as those used 
in ICooper et al.l (|2006f) . ther efore our Figu r e [T^ d oes not 
look identical to Figure 5 in ICooper et ah! (|2006[ ). 

The overall agreement is very good between the rela- 
tive overdensity and the relative bias from the cluster- 
ing measures presented here; these two statistics should 
be similar but not identical. The overdensity estimate, 
(^3, is not measured at a fixed scale but is generally in 
the range r = 1 — 3 h^^ Mpc; we compare to the rela- 
tive bias from clustering averaged on scales = 1 — 5 
Mpc to minimize systematics from scale dependence. 
The environment overdensity statistic is in some ways a 
scale-averaged clustering measure, so the agreement be- 
tween the two should be good. However, the environment 
statistic uses all neighboring galaxies to define the local 
overdensity, not just those with the same color or magni- 
tude as the galaxy for which one is measuring the envi- 
ronment. This is different from the clustering statistic we 
use, where only galaxies within a given color and mag- 
nitude range are used for the auto-correlation function. 
Using all galaxies as a density tracer could potentially 
cause the environment statistic to be more sensitive to 
group membership and corresponding small-scale over- 
densities for populations which tend to be central galax- 
ies instead of satellite galaxies (i.e., if there is only one 
galaxy of that type in an overdensity). 

As seen in Figure [121 the agreement between the two 
statistics is good but not exact. In particular, the clus- 
tering results show a continued trend within the blue 
population of the bluest galaxies being the least clus- 
tered, which is not clearly reflected in the environment 
measure, though the discrepanc y is at the la level; the 
broader environment samples in ICooper et al.l (|2006f ) do 
show such a downturn. 

Our clustering-based results on the color-density rela- 
tion at z 1 are very simi l ar to l ocal e nvironment find- 
ings in SDSS. iHogg et all (|20M l200l find that SDSS 
galaxies also show a color-density relation within the blue 
population, but not within the red sequence. They also 
find that locally the overdensity of blue galaxies shows 
little luminosity dependence, while there is strong lu- 
minosity dependence within the red galaxy population. 
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Unfortunately, the DEEP2 volume is not large enough to 
provide a fair sample of the brightest, rarest red galaxies 
at z ~ 1, and the error bars on our results for tq as a 
function of luminosity within the red sequence are 5-15%, 
which may be large enough to mask a modest luminosity 
dependence for the range of luminosities we probe. The 
general conclusion is that there has not been qualitative 
evolution in the color-luminosity-overdensity space since 
z - 1. 

4.8. Pairwise Velocity Dispersion 

On small scales {r ^ I h^^ Mpc) random motions 
within galaxy groups and clusters lead to elongations of 
the ^(rp, tt) contours along the line of sight, the so-called 
"fingers of god" . These motions can only be seen in spec- 
troscopic redshift surveys, and at z ~ 1 only the DEEP2 
and VVDS surveys can currently detect these distortions 
in ^(rp,7r). The virialized motion of the galaxies can 
be measured using the pairwise velocity dispersion, ai2, 
which is estimated by modeling ^(r^, tt) as a convolution 
of the real-space correlation function, ^(r), with a broad- 
ening velocity function. The relation between £,{rp,Tr) 
and ^(r) is usually assumed to be 

l + e(rp,7r)= J f{vi2) [l + ^i^r^+y{vi2,n)^)] dvu, 

(14) 

where /(W12) is the distribution function of the relative 
velocity dif ference between galaxy pairs along the line of 
sight (e.g.. iDavis fc Peebles! 11983 ). This velocity differ- 
ence is defined as 

V12 = TT - y + vi2{r) (15) 

where V12 (r) is the mean radial pairwise velocity of galax- 
ies at separation r and y is the real-space separation 
along the line of sight. Here y has units of km s~^ in- 
stead of Mpc (this conversion includes a factor of 
Hq, with h = 1, and a factor of (1 + z)), such that the 
same physical velocity aX z — \ has less of an effect on 
£,(Tp, tt) than at z = 0.1. 

To measure the small-scale velocity disp e rsion o n scales 
r ~ 1 Mpc, we follow iFisher et all (| 19941) and in 
Equation [13] we use 

^(rp = l,7r) = 0.5[C(rp = O.SMpc /i'^ 7r)+^(rp = 1.5Mpc/i^ 

(16) 

in both real and redshift space, for values of tt < 20 
Mpc. We then normalize ^(rp = l,7r) so that our 
subsequent fitting will be sensitive to the overall shape 
of ^(rp = 1, tt) but insensitive to the amplitude. 

An e xponent ial form is usually adopted for fjvjj) 
/e.g.. IDavis fc Peebles 1983: [Fisher et all [l994 : 
iDiaferio fc Gerieiiil996i: iShetbil99l : 

/(«12) = CXp [ \vi2 - iJl2| I , 

V2(T12 V Cri2 I 



to the assumed model; (T12 is ^--^20-40% lower if TJ12 is 
assumed to be negligible. 

To estimate CT12, £,{r) is modelled using the power-law 
fits for To and 7 for each sample, and we minimize 
between the observed ^(rp = l,7r) and modelled ^(rp — 
l,7r) for a range of (T12 values. The results are given 
Table 2, where errors are derived from jacknife samples. 
Red galaxies are found to have au — 530 ± 50 km s~^, 
while blue galaxies have a significantly lower value of 
(J12 = 240 ± 20 km s~^. The error bars are too large to 
detect significant differences as a function of luminosity 
within the red or blue galaxy populations or as a function 
of color for the finer color bins. Green galaxies have a 
large value, ai2 — 490 ± 110 km s~^, consistent with red 
galaxies and 2a larger than blue galaxies. 

These numbers are similar to, but somewhat smaller 
than, values found at z = 0.1. Using the same infall 
model. IZehavi et all (|2002f) measure ai2 ~ 650 — 750 km 
s~^ for red galaxies in SDSS and ai2 ^ 300 — 45 km 
s~^ for blue galaxies at z = 0.1. iLi et al.l (|2006b( l also 
find for SDSS galaxies that ai2 ~ 600 - 800 km s'^ 
for red g alaxies and ai2 ^ 200 — 400 km s~^ for blue 
galaxies. iMadgwick et al.l ()2003af ) measure similar values 
for galaxies in the 2dF Redshift survey: ai2 = 612 ± 92 
km s^^ for passive galaxies (defined by spectral type) and 
<7i2 — 416 ± 76 km for active, star-forming galaxies. 

Our finding that red galaxies have stronger fingers of 
god at z ~ 1 than blue galaxies again reflects the result 
that, on average, red galaxies reside in more overdense 
regions, such as galaxy groups, than blue galaxies. 

4.9. Multipole Moments 

Peculiar velocities also affect C(rp,7r) on large scales, 
where coherent infall of galaxies onto forming structures 
flattens ^(j'pjTr). A standard method for quantifying 
these large-scale redshift-space distortions is to measure 
the ratio of the quadrupole to r nonopole mome nts of 
the two-point correlation function (lH amiltolilll992| ). The 
two-dimensional correlation function can be decomposed 
into a sum of Legendre polynomials, 



i{rp,7r) = }_^Ciis)Viifi), (19) 

\-)] 

where Vi is the Legendre polynomial and /i is the co- 
sine of the angle between the line of sight and the redshift 
separation vector, s. The multipole moments are defined 
as 



6(s) = 



2/ + 1 



(17) 



Following iHamiltmil ()1992f ). we define 

6(5) 



Q{s) 



(20) 



(21) 



where V12 is the mean and (T12 is the dispersion of the 
pairwise peculiar velocities. We assume an infall model 
based on the similarity solution of the pair conservation 
equation ([Davis fc Peeble£l[T977l) . 

where r^, = 5 h^^ Mpc and y is the radial separation 
in the real space. The results are relatively insensitive 



This ratio depends on / 3 = /b, where h is the linear 
galaxy bias ([Kaiseij|1987D . as 



Q{s) = - 



(4/3)/3+(4/7)/3^ 
l-t-(2/3)/3-H(l/5)/?2 



(22) 



Figure [Til shows Q{s) for the main blue and red galaxy 
samples. Q{s) is positive on small scales where the fin- 
gers of god are strong and negative on large scales where 
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coherent infall dominates. The dotted lines indicate Q 
predicted on linear scales (from Equation 22), assuming 
Vim at the mean redshift of each sample (for the blue 
sample 51™ (z = 0.90) = 0.75 and for the red sample 
^ra{.z = 0.82) = 0.72, corresponding to rim(z = 0) = 
0.3) and using the inferred linear galaxy bias for each 
sample, h = 1.28 for blue galaxies and b = 1.65 for red 
galaxies. Using flm{z = 0) = 0.24 increases the predicted 
Q values by ~5%. For the blue sample (3 = 0.66 and for 
the red sample (3 = 0.50. Both galaxy samples have neg- 
ative Q values, on scales s > 5 Mpc for blue galaxies 
and s > 10 Mpc for red galaxies, indicating coher- 
ent infall, though the significance is lower for red galax- 
ies. On the largest scales Q increases; however, the error 
bars increase on these scales as well, which approach the 
transverse scale of the survey and so may be susceptible 
to systematic effects. On scales s ~ 10 — 15 Mpc the 
blue sample is consistent with the predicted constant Q 
value within the noise. The red sample does not reach the 
predicted negative value on the largest scales; however, 
given the non-negligible pairwise velocity dispersion we 
do not expect this model to be perfect, as even on scales 
of 15 Mpc, Q{s) is sensitive to ai2 and the form of 
fivu)- This is especially true for the red galaxy sample, 
where ai2 is greater. The fact that infall is seen indicates 
that the clustering of galaxies is dynamic and growing, 
as expected. 

5. COMPARISON WITH SEMI-ANALYTIC MILLENIUM 
RUN SIMULATION 

We compare our clustering results for blue and red 
galaxies with the recent sem i-analytic galaxy formation 
model of ICroton et al.l ()2006[ ) applied to t he Millennium 
Run N-body dark matter simulation (jSpringel et al.l 
l2005f ). The Millennium Run follows the dynamical evo- 
lution of 10^" dark matter particles in a periodic box of 
side- length 500 h^^Mpc with a mass resolution per parti- 
cle of 8.6 X 10* h~^MQ. Cosmological parameters consis- 
tent with those measured from the first year WMAP data 
(|Spergel et al.l[200l ISelTak et all 12005 1) were used. The 
galaxy formation model follows the growth of approxi- 
mately 25 million galaxies from their birth to the present 
day and has been tuned to provide a good match to many 
observed properties of local galaxies, primarily the lumi- 
nosity function and TuUy-Fisher relation. Within this 
model the behavior of the observed global galaxy two- 
point correlation function is well reproduced at z = 0, 
however when split by galaxy color, z = red galax- 
ies are somewhat too cluster ed in the model whi le blue 
galaxies are under-clustered (Spring el et al.|[2005f) . Ours 
is the first test of the clustering of Millennium Run galax- 
ies as a function of color at higher redshift. 

Using an ensemble of 15 light-cones, each of which 
matches the geometry of one of the fou r DEEP2 sur- 
vey fields (see iKitzbichler fc White! I2007L for further de- 
tails) we measure the clustering of blue and red galaxies 
with the same Mb and z ranges as our main DEEP2 
color samples. We define galaxies to be blue or red in 
the model using the same color cut as used for DEEP 2 
galaxies (Equation 1). If we instead use the minimum 
of the observed color bimodality in the model galaxies 
(which is roughly 0.1 mags blueward of the location for 
the DEEP2 galaxies) our qualitative results and conclu- 
sions remain unchanged. Details of each sample used in 



the mock catalogs are given in Table 2. The number den- 
sities for red and blue galaxies in the model match the 
observed number densities quite well. 

The resulting ^(rp,7r) contours for the semi- analytic 
model galaxies are shown in Figure[T5l Red model galax- 
ies are more clustered and have stronger fingers of god 
than the blue model galaxies. A direct comparison be- 
tween the red model galaxy ^(7'p,7r) and the equivalent 
DEEP2 result (Figure [3|) indicates that a) the redshift- 
space distortions are generally similar between the model 
and data, including both the fingers of god and the co- 
herent infall on large scales, and that b) the overall clus- 
tering amplitude is larger in the model; e.g. compare the 
thick ^ = 1 lines. Hence, the discrepancy seen at z = 
appears to extend back to earlier epochs. The blue model 
galaxy S.irp,!:) shows differences with the data as well. 
The left panel of Figure [15] shows coherent infall on large 
scales but no small-scale redshift-space distortions, un- 
like the DEEP2 data. This is reflected in the low value 
of (Ti2 found for the model blue galaxies (see Table 2). 

The discrepancies seen in Figure [T5] are highly sig- 
nificant, as seen in the projected correlation functions 
shown in Figure \W[ where the observed DEEP2 results 
are shown for comparison. On scales rp < 0.5 Mpc 
the blue cloud semi-analytic model correlation function 
is flat and does not continue to rise as in the data. In ad- 
dition, in the quasi- linear regime, rp — 1 — 10 Mpc, 
the blue model galaxies are less clustered than the data 
by ^ 15%. In contrast, the red galaxies in the model are 
significantly more clustered than the data, ~ 40% over 
most scales, as also seen at low redshift. Power-law fits 
are given in Table 2, where slightly larger scales are used 
for the model fits {rp = 0.4 - 20 Mpc) than for the 
data, as the model blue galaxy Wp{rp) is not a power law 
on smaller scales. 

Finally, the red-blue cross-correlation function for the 
model galaxies is shown in Figure [17] (solid line) along 
with the cross-correlation function in the DEEP2 data 
(thin dashed line) and the geometric mean of the red 
and blue auto-correlation functions in the model (dot- 
dash line). The cross-correlation function in the model 
is clearly lower on small scales (r^ < 1 h^^ Mpc) than 
either the geometric mean or the DEEP2 data. This 
refiects the lack of blue galaxies around red galaxies in 
the model. Interestingly, a hint of this same fiattening 
of Wp(rp) for blue g alaxies in this model can be seen at 
z = in Figure 5 of ISpringel et al.l (|2005l ): however they 
do not plot below scales rp = 0.5 Mpc where for us 
the effect is most pronounced. 

It appears that the discrepancies shown in Figures \W\ 
and [17] are due to the incorrect modeling of the physics 
that determines the colors of satellite galaxies. If too 
few model satellite galaxies are blue compared to the 
data at z '--^ 1, then both the model blue auto-correlation 
function and red-blue cross-correlation function will have 
lower amplitudes on small scales than seen in the DEEP2 
data, inside the 'one- halo' regime {rp < 1 h^^ Mpc). 
Similarly, an overabundance of red satellites on such 
scales should also be present, as is seen in Figure [TBI 

A similar effect is detected at low redshift when 
comparing the fraction of blue SDSS galaxies in 
groups with thes e Millenium Run semi- analytic models. 
IWeinmann et al.l (|2006a ) find that the blue fraction of 
galaxies in groups is much higher in the SDSS data than 
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in the model, which predicts at z = 0.1 that almost all 
satellite galaxies (~85%) are red. 

Our results provide a valuable way to test some of the 
physical assumptions made by the semi-analytic model 
to better understand the processes that may actually 
be occurring and their relative importance. We iden- 
tify two overly simplified aspects of the model satellite 
evolution that contribute to the clustering discrepancies 
shown in Figure [TBI The first is the way in which hot gas 
is stripped from halos (Mhaio ^ IO^^H'^Mq) when the 
halo falls into t he potential of a more m assive system. 
As discussed bv IWeinmann et all (|2006al ). the common 
assumption in most semi-analytic models is one of ex- 
treme strangulation, where all of the hot gas is instantly 
stripped and added to that of the larger halo upon ac- 
cretion. This simplifying assumption is important for 
galaxy colors because the dynamical time for a sub-halo 
(and hence satellite galaxy) to merge is usually longer 
than the time it takes the stellar population to fade and 
redden, which all model satellites will do in the absence 
of fresh fuel for star formation. This provides an overly- 
efficient channel for blue satellites to transform into red 
s atellites. 

IWeinmann et all ()2006aD point out that the 
iCroton et a l.' ('2006') model does not include a host 
halo mass dependence for the strangulation efficiency. 
Given that the ability of a parent halo to strip gas 
from a subhalo likely depends on the density of the 
hot gas in the parent halo, a model which includes 
such a mass-dependent strangulation prescription might 
provide a better match to the results shown here. 

The second aspect of the model satellite evolution that 
may be discrepant with the real universe is the assump- 
tion that local potential wells (i.e. subhalos inside larger 
parent halos) cannot accrete gas, meaning that after a 
satellite galaxy is captured by the larger halo it receives 
no new baryonic fuel for later star formation from the 
intragroup or intracluster medium. This may be gener- 
ally true in N-body simulations which show that subhalos 
tend to lose mass through tidal stripping as they spiral 
inward, not gain mass. However it is plausible that in the 
outer regions of a dark matter halo, the local potential 
of a subhalo may dominate that of its host and attract 
new baryons from the hot parent halo, baryons that may 
subsequently cool onto the satellite. Such gas would ex- 
tend the star formation history of the satellite and hence 
modify the evoution of its color. 

It is also important to note that there are processes 
that are known to operate in cluster (and also group) en- 
vironments that act to remove cold gas from a satellite 
galaxy, such as tidal stripping an d harassment. These 
processes are not included in the ICroton et al.l ()2006f ) 
semi-analytic model used here (most semi-analytic mod- 
els ignore such additional effects to keep the models rel- 
atively simple). Their inclusion would typically redden 
the color of the satellite further with time, worsening an 
effect that is already too efficient in the current model. 
Our work indicates that current prescriptions for satellite 
infall and its consequent effect on star formation may be 
too simplistic. However, it also highlights the utility of 
color-dependent clustering measurements in constraining 
such models. 

6. SUMMARY AND DISCUSSION 



In this paper we use volume-limited subsamples of the 
DEEP2 Galaxy Redshift Survey data to measure the 
color and luminosity dependence of galaxy clustering at 
z ^ 1. We split the sample into blue, star-forming galax- 
ies and red, nonstar-forming galaxies using the observed 
color bimodality in {U — B). We further subdivide each 
sample into luminosity and finer color bins. 

In this section we first summarize our main results and 
then discuss their implications. Our main findings are: 

• Red galaxies are much more strongly clustered than 
blue galaxies at z '--^ 1, with a relative bias that is 
roughly as high as is found locally {bred/bbiue ~ 
1.4 — 1.6). This implies that, for relatively bright 
galaxies, the color-density relation was as strongly 
in place 8 Gyr ago as it is observed to be today. 

• The comoving clustering amplitudes that we mea- 
sure for both blue and red galaxies are only 
slightly (10-20%) lower than what is found lo- 
cally. The measured values are somewhat higher 
than those in prev i ous DEEP2 and VVD S results 
(|Coil et al.l l2004bl : iMeneux et all I2006D . though 
generally within 2a. 

• We find no statistically-significant dependence of 
galaxy clustering on luminosity for either red or 
blue galaxies, within the luminosity range that 
we probe. Additionally, within the red sequence 
there is no dependence of clustering on color, while 
within the blue cloud there is a strong dependence 
on color. The color dependence and lack of lumi- 
nosity dependence for both red and blue that we 
find is consistent with previous results at z = 0.1, 
for the luminosity range probed here. 

• The stronger dependence of clustering on luminos- 
i ty found for all galaxies in the DEEP2 sample 
()Coil et al.l2006af ) is likely due in part to the chang- 
ing fraction of red versus blue galaxies as a function 
of luminosity, as we do not find as strong of a lu- 
minosity dependence here within either the red or 
blue galaxy populations independently. 

• The brightest blue galaxies have a significant rise in 
their correlation function on small scales (j-p < 0.2 

Mpc), which likely reflects a lower satellite 
fraction compared to fainter blue galaxies; i.e., 
bright blue galaxies are more likely to be central 
galaxies within their parent dark matter halos. 

• The correlation between color and clustering within 
the blue cloud is likely a reflection of a correlation 
between stellar mass and halo mass. The redder 
galaxies in the bl ue cloud have higher s tellar masses 
(see Figure 8 of ICooper et al1l2007bD and are sig- 
nificantly more clustered than the bluer galaxies, 
which implies that they are in more massive dark 
matter halos. 

• Red galaxies show strong 'fingers of god' at z ~ 1, 
indicating that they lie in virialized overdensities 
such as groups and clusters. Blue galaxies show 
smaller 'fingers of god'. 
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• Both blue and red galaxies show a flattening of 
^(rp,7r) on large scales due to coherent infall of 
galaxies onto structures that are still collapsing 
(the Kaiser effect). 

• The projected correlation functions of brighter 
samples near L* show deviations from a power law 
on small scales, within the one-halo regime. The 
deviations are more pronounced for blue galaxies 
than for red galaxies. 

• The red-blue cross-correlation function is con- 
sistent with the geometric mean of the auto- 
correlation functions of blue and red galaxies sep- 
arately. There is a marginal trend on small scales 
for the cross-correlation to lie below the geometric 
mean. Such a feature, if real, would indicate a lack 
of blue galaxies near red galaxies. This could reflect 
a suppression of star formation in the centers of 
overdensities, such as groups, which are dominated 
by red galaxies at z ~ 1 as seen in the group- galaxy 
cross-correlation function ijCoil et al.ll2006lj) . 

• Green galaxies in the valley of the observed color 
bimodality display a similar large-scale clustering 
strength as red galaxies, but show a small-scale 
clustering amplitude and infall kinematics akin to 
blue galaxies. The green galaxy population does 
not have strong 'fingers of god' but does show co- 
herent infall on large scales, and the correlation 
function has a similar slope to the blue galaxy pop- 
ulation. Green galaxies thus appear to be at the 
edges, but not in the centers, of the same over- 
densities as red galaxies. This is consistent with 
green galaxies being transition objects moving from 
the blue cloud to the red sequence and supports a 
picture of quenching happening preferentially on 
the outskirts of overdensities at z '--^ 1 on scales 
Tp -0.5-2 /i-i Mpc. 

• We provide a direct test of environment and clus- 
tering statistics, comparing the relative bias as a 
function of color using a local galaxy overdensity 
estimator with the two-point correlation function 
results, and find good agreement between the two. 

• Comparing wi th the semi-analytic galaxy evolu- 
tion model of iCroton et al.l ()2006f ). we find that 
red model galaxies are significantly more clustered 
than the DEEP2 galaxies, and that in the model 
the clustering strength of blue galaxies is too low, 
especially on small scales {rp < 0.5 h^^ Mpc). 
These discrepancies are likely due to galaxies be- 
ing quenched too efficiently in the model; i.e., star 
formation is shut down in satellite galaxies too 
quickly. 

• These results all rely on precise, robust, spectro- 
scopic redshifts. The DEEP2 rms redshift errors 
(determined from repeated observations) are < 35 
km s~^, which is unprecedented at z ~ 1 and makes 
the ^(^p, tt) analyses presented here possible. 

From the observed bias of the main red and blue 
galaxy samples we can infer the minimum dark mat- 
ter mass for halos that host these galaxies. Un- 



der the assumption that all galaxies are central galax- 
ies and not satellites (i.e., one galaxy per dark mat- 
ter halo), the large scale bias of the galaxy popula- 
tion can be matched to the bias of dark matter ha- 
los using the formulae from e.g., Mo & White (1996) 
or Sheth, Mo & Tormen (2001). However, we know 
that not all of the DEEP2 galaxies are central galaxies 
(iGerke et a l."2005': 'Co ^ov.JVechsler. fc Kravtsovl [20061 : 
iZheng. Coi l. & Zchavi 20071 ). and to obtain a more 
precise estimate of the minimum dark matter halo 
m ass we use the results pres ented in the Appendix 
of IZheng. Coil, fc Zehavil (|2007). which take into ac- 
count the inferred satellite fraction from HOD model- 
ing of the DEEP2 lumi nosity-dependent clustering re- 
suhs ()Coil et all |2006a'). At z = 0.9, the observed 
bias and number densities of the main red and blue 
galaxy samples (with Mb < —20.0) imply minimum 
dark matter halo masses of Mmin ~ 2 x IQ^'^h^^MQ 
for red galaxies and Mmin ~ 5 x 10^^/i~^Mq for 
blue galaxies. However, the halo occupation of red 
and blue galaxies is likely complicated enough to war- 
rant a full HOD model, which we defer to a fu- 
ture paper (see, e.g.. Ivan den Bosch, Yang, fc Mo|[2003l : 
ICollister fc Lahavll2005l : IZehavi et al.ll2005l for ^dehng 
of local red and blue samples). 

The different slopes measured for the clustering of red 
and blue galaxies would pose a potential difficulty in 
terms of understanding how red and blue galaxies popu- 
late dark matter halos, if the slopes were different on 
large scales (r^ > 2 Mpc). On these scales the 
halo bias is linear, in that there is an amplitude offset 
between ^(r) for halo s of different mass, but no differ- 
ence in the slope (e.g.,lvaji den Bosc h. Yang, fc Mo|[2003l : 
iConrov. Wechsler. fc Kravtso v 2006). We do not have a 
large enough survey here to claim a significant difference 
in the slope on scales > 2 Mpc, though the upper 
left panel of Figure H] does suggest such a difference. 

It has been previously established that the num- 
ber density of gala xies on the red sequence has in- 
creased since z = 1 (jBefl et al.l[200l iFaber et al.ll2007l : 
iBrown et al.l l2007h . It is thought that these galaxies 
likely migrated from the blue cloud to the red sequence 
as their star formation was quenched. Using the observed 
clustering of red and blue galaxies at z ~ 1 , one can pre- 
dict what their clustering strength should be at z ~ 0, for 
a given cosmology and assuming no merging of galaxies, 
by using the contin uity equation and the exp ected growth 
of perturbations (Tegmark fc Peebieill998f ). This leads 
to larger predicted clustering amphtudes than is observed 
for local L* red and blue galaxies; the main red and blue 
samples used here (both of which have L ~ L* at z = 0.9) 
would have clustering scale- lengths of ro = 6.5 Mpc 
and ro = 5.5 h^^ Mpc, respectively, at z = 0. The fact 
that local clustering measurements at L* are lower than 
these predictions, especially for blue galaxies, lends cre- 
dence to the idea that blue galaxies have migrated to the 
red sequence and less clustered blue galaxies have joined 
the blue cloud since z = 1. 

The color and luminosity dependence of clustering at 
z = 1 can illuminate which blue galaxies likely become 
red between z = 1 and z = 0. The fact that the relative 
bias between red and blue galaxies is as high at z ~ 1 
as z '--^ indicates that blue galaxies have to be turning 
red in a density- dependent way, in that the more clus- 
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tered of the blue galaxies must become red. The bias for 
any R alaxy pop ulation should tend towards unity with 
time (jFrvlllQQa ) and linear growth theory would predict 
that the relative bias of our red and blue samples should 
decrease by 8% to z = 0.1. Observationally, however, 
it appears that the relative bias between red and blue 
galaxies has either remained constant since z ^ 1 or in- 
creased. This therefore implies that the most clustered 
blue galaxies must be turning red to boost the clustering 
of red galaxies relative to blue galaxies and to counteract 
the effect of the relative bias decreasing. It is the most 
clustered of the blue galaxies at z = 1 that are likely 
to be as clustered as red galaxies by z = 0, and our re- 
sults here indicate that the reddest of the blue galaxies 
and the green galaxies are the most likely candidates for 
having their star formation quenched by z = 0, based on 
their clustering properties. 

The luminosity dependence of clustering at z ~ 1 can 
be explained as being almost entirely due to the rela- 
tionship between luminosity an d halo mass, as shown by 
IConrov. Wechsler. fc Kravtsovl p006) who use a simple 
one-to-one relation between halo mass and light to obtain 
the same luminosity- and scale-dependence of clustering 
in an N-body simulation as we find in DEEP2. The tight 
relation between halo mass and galaxy light is confirmed 
by Zheng, Coil. & Zehavi (2007), who fit HOD models 
to the DEEP2 luminosity dependent clustering results. 
However, the color dependence of clustering seen here 
can not entirely be halo-mass dependent, as the cluster- 
ing as a function of color does not match the results of 
IConrov. Wechsler. fc Kravtsovl ([2006), where the galax- 
ies have been ranked by luminosity, which is equivalent 
to halo mass; i.e., the trends seen with color in the data 
do not match the trends found with luminosity or halo 
mass in the model. This implies that there is a wider 
range in the Mhaio/L ratio of galaxies at a given color 
than at a given luminosity. To constrain galaxy evo- 
lution physics it may therefore be more informative to 
study samples divided by color rather than luminosity, 
as more galaxy physics is apparently needed to interpret 
the color dependence of galaxy clustering. The lack of 
a stronger luminosity dependence for the separate red 
an d blue galaxy samples shown here ma kes the results 
of IConrov. Wechsler. fc Kravtsovl (|2006f ) more surpris- 
ing, as it is not clear why Mb luminosity for all galaxies 
should be so tightly correlated with halo mass. It may be 
that the relative mix of red and blue galaxies as a func- 
tion of luminosity helps to create a stronger relation with 
luminosity for all galaxies and is more tightly correlated 
with halo mass than the luminosity itself. 

As noted in previous papers o n the color-density 
relation in the p EEP2 sample ()Gerke et all 120071 : 
ICooper et al.ll2006l ). the observed color-density relation 
at z ~ 1 can not be the result of cluster-specific physics 
as only a few percent of DEEP2 galaxies are in massive 
clusters. The color dependence of clustering must be 
dominated by either processes that operate in a group 
environment or intrinsic processes that depend on the 
age, stellar mass or halo mass of a galaxy. Processes that 
occur when a galaxy is accreted into a larger halo, be- 
coming a satellite galaxy, may contribute to the observed 
color-density relation but do not dominate it, as most 
DEEP2 galaxies are not satellite galaxies (jCerke et al.l 
l2005HZheng. Coil, fc Zehavill2007f l. What star-formation 



quenching processes may be at work in groups? There 
are more galaxy mergers in groups than in clusters, such 
that mergers may play an important role in shutting off 
star formation in group galaxies. There is less hot ambi- 
ent gas in groups than in clusters; it is not clear whether 
the density of the intragroup gas is high enough to strip 
gas from infalling satellite galaxies, except near the core 
of the group. Apparently not all satellite galaxies are red, 
as is seen in the clustering results; this implies that there 
is a non-negligible timescale for quenching star formation 
in satellite galaxies. 

DEEP2 galaxies in groups at z '--^ 1 have a clustering 
scale-length of r n = 4.97 ± 0.25 Mpc and a slope of 
7 = 2.15 ± 0.06 (|Coil et al.ll2006bD . similar to red galax- 
ies as measured here. Galaxies in the field at z '--^ 1 have 
a slope of 7 = 1.70 ± 0.13, similar to blue galaxies, but 
with a lower correlation length, rg — 2.54 ± 0.25 
Mpc. The cross-correlation of red galaxies an d groups is 
higher than that of blue galaxies and groups (jCoil et al.l 
l2006b( l. which also indicates that red galaxies are more 
likely to be in groups. However, the picture is not as 
simple as all red galaxies being in groups and blue galax- 
ies being in the field: blue galaxies have some fingers 
of god on small scales, indicating that at least some are 
in virialized overdensities, and the brighter blue galaxies 
show a rise in the correlation function on small scales, 
which is probably due to a lower satellite fraction for the 
brighter blue samples. These galaxies are likely the cen- 
tral galaxies in groups at z = 1, and may migrate to the 
red sequence by z 0. The scale at which the one-halo 
term rise occurs for these samples is small, Vp ^ 0.3 
Mpc (as seen in Figure HO]) , and indicates again that these 
galaxies are likely to be in the centers of the groups they 
reside in. Additionally, ~ 5 — 10% of DEEP2 gal axies in 
underdense regions are red (|Cooper et al.ll2007aD . 

The green galaxy population, which lies near the min- 
imum in the observed color bimodality, appears to be, at 
least in part, a distinct population and is not a simple 
mix of red and blue galaxies. Though the green popula- 
tion can include the tail of each of the blue and red distri- 
butions, it also contains transitional objects, likely mov- 
ing from the blue cloud to the red sequence. The green 
galaxies have different clustering characteristics than ei- 
ther the red or blue galaxies as a whole. On large scales 
[vp > 1 Mpc) the green galaxies are as clustered as 
red galaxies of the same luminosity, implying that they 
reside in the same relatively massive halos as red galax- 
ies. On small scales, however, the clustering amplitude 
is similar to that of blue galaxies, which likely reflects 
a lower radial concentration for green galaxies compared 
to red galaxies. The general picture that emerges is one 
in which star formation quenching is occurring primarily 
at the edges of overdensities as galaxies fall in, turning 
blue galaxies green around the edges of groups. 

We have identified that the scale at which this tran- 
sition from red to blue occurs at z ~ 1 is 0.5 h^^ Mpc 
< Tp < 2 Mpc. This scale is just larger than the tran- 
sition scale between the o ne-halo and two-halo terms at 
z ~ 1, which is shown bv IZheng. Coil, fc Zehavi! ()2007D 
(see th eir Fig. 1) to be ^ 0.5 Mpc for all galaxies 
and bv lCoil et all ()2006bl ) (see their Fig. 8) to be - 1.0 
Mpc for galaxies in groups. This transition scale for 
green galaxies is therefore similar to or just larger than 
the largest halos, which indicates that it is occurring at 
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~ 1 — 2 virial radii, on the edges of large halos. This 
conchision lends support to the idea that the mechanism 
which shuts off star formation in these galaxies is not 
happening at the centers of massive halos, but on the 
outskirts. 

The comparison to semi-analytic galaxy formation 
models demonstrates vividly the utility of using corre- 
lation statistics to test theoretical models, and the need 
to test them at intermediate redshift, not just locally. 
It is also likely that currently popular models that have 
a simple halo threshold mass for quenching star forma- 
tion in galaxies will suffer similar problems as the model 
shown here. We urge theorists to use these data to con- 
strain their models. 
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TABLE 1 

Published Color Clustering Results 





mean 


blue 


blue 


red 


red 


relative 


M 






z 




7 




7 


bias 


range 




Norberg et al. 2002 


0.1 


4.45 ± 0.47 


1.76 ± 0.09 


5.71 ±0.57 


1.87 ±0.09 




-20 < M^j < 


-19 


Madgwick ct al. 2003 


0.1 


3.67 ±0.30 


1.60 ± 0.04 


6.10 ±0.34 


1.95 ±0.03 


1.45 ±0.14 






Budavari ct al. 2003 


0.2 


~4.5-6.0 


~ 1.67 ±0.10 


~ 6.3 - 6.6 ±0.20 


2.0 ±0.5 




-20 > Mr, > 


-21 


Zchavi ct al. 2004 


0.1 


3.63 ± 0.16 


1.69 ± 0.04 


5.67 ±0.37 


2.08 ±0.05 


~ 1.6 


-20 < Mr < - 


-19 


Coil ct al. 2004 


0.8 


2.81 ± 0.48 


1.52 ± 0.06 


4.32 ±0.73 


1.84 ±0.07 


1.41 ±0.10 






Meneux et al. 2006 


0.8 


2.18 ±0.30 


1.40 ±0.14 


3.78 ±0.7 


1.87 ±0.25 


1.45 ±0.27 







TABLE 2 
DEEP2 Galaxy Samples 



Sample 


No. of 
galaxies 


n 

{h^Mpc~ 




median 
(U-B) 


Mb 
range 


Median 
Mb 




z 

range 




Mean 
z 








7 




Bias"* 




cri2 


Red: -21.0 


792 


3.9 


10- 


-4 


1.22 


< 


-21.0 


-21.38 





7 


- 1 


1 


0.92 


5 


78 ± 0.87 


1 


89 ±0.17 


1 


79 ±0 


21 


360 ± 70 


Red: -20.5 


1370 


9.0 


10" 


-4 


1.21 


< 


-20.5 


-21.03 





7 


- 1 


025 


0.88 


5 


37 ±0.27 


1 


94 ± 0.06 


1 


71 ±0 


14 


500 ± 50 


Red: -20.0 (main) 


1474 


1.6 


10 


-3 


1.21 


< 


-20.0 


-20.70 





7 


- 


925 


0.82 


5 


25 ± 0.26 


2 


06 ± 0.04 


1 


65 ±0 


15 


530 ± 50 


Red: -19.5 


1032 


2.4 


10 


-3 


1.20 


< 


-19.5 


-20.46 





7 


- 


825 


0.77 


5 


04 ± 0.40 


2 


06 ±0.12 


1 


55 ±0 


15 


490 ± 40 


Blue: -21.0 


1424 


6.1 


10 


-4 


0.81 


< 


-21.0 


-21.28 





7 


- 1 


2 


1.00 


4 


27 ±0.43 


1 


75 ± 0.13 


1 


39 ±0 


16 


250 ± 50 


Blue: -20.5 


3747 


1.6 


10 


-3 


0.77 


< 


-20.5 


-20.89 





7 


- 1 


2 


0.99 


3 


86 ± 0.15 


1 


77 ±0.05 


1 


30 ±0 


07 


360 ± 30 


Blue: -20.0 (main) 


4808 


3.5 


10 


-3 


0.74 


< 


-20.0 


-20.48 





7 


- 1 


05 


0.90 


3 


87 ±0.12 


1 


64 ± 0.05 


1 


28 ±0 


04 


240 ± 20 


Blue: -19.5 


4018 


6.3 


10 


-3 


0.68 


< 


-19.5 


-20.10 





7 


- 


9 


0.81 


3 


58 ± 0.25 


1 


66 ±0.07 


1 


24 ±0 


08 


270 ± 20 


Red: Redder 


731 


7.9 


10 


-4 


1.27 


< 


-20.0 


-20.66 





7 


- 


925 


0.81 


5 


12 ± 0.51 


2 


17 ±0.07 


1 


60 ±0 


21 


740 ± 180 


Red: Bluer 


743 


8.1 


10 


-4 


1.12 


< 


-20.0 


-20.61 





7 


- 


925 


0.82 


5 


18 ± 0.57 


2 


13 ± 0.09 


1 


61 ± 


22 


740 ± 190 


Green 


997 


7.7 


10 


-4 


1.00 


< 


-20.0 


-20.65 





7 


- 1 





0.87 


5 


17 ±0.42 


1 


59 ± 0.08 


1 


65 ± 


26 


490 ± 110 


Blue: Redder 


1407 


1.3 


10 


-3 


0.88 
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-20.0 


-20.65 





7 


- 1 





0.88 
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49 ± 0.35 
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64 ±0.14 
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44 ±0 


19 


320 ± 40 


Blue: Middle 
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1.1 
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0.72 
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-20.0 


-20.44 
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- 1 
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63 ± 0.18 
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84 ±0.12 
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270 ± 50 


Blue: Bluer 
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1.1 
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0.56 
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-20.33 
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- 1 
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16 ± 0.15 
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78 ±0.14 
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11±0 
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260 ± 70 


MR+SAM Red^: 
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-3 


1.14 
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-20.0 


-20.40 





7 


- 1 





0.89 
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31 ± 0.25 


1 


91 ±0.04 


1 


96 ±0 


14 


490 ± 30 


MR+SAM Blue: 


23645 


4.2 


10 


-3 


0.65 
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-20.0 


-20.43 





7 


- 1 





0.88 


3 


27 ± 0.07 


1 


50 ± 0.03 


1 


20 ±0 


10 


140 ± 20 



^ For a ACDM model with cts = 0.9. 
MR+SAM: Millenium Run with semi-analytic model 
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Fk;. 1. — Top row: Rcstfrarnc color-magnitude diagrams for DEEP2 galaxies in the redshift range 0.7 < z < 1.0. The bold vertical 
dashed line indicates the Mg = —20.0 luminosity cut used for our main color samples, while the tilted and horizontal dashed lines show 
the various color divisions used to create color subsamples. The left panel shows the cuts for the main red and blue samples (bold dashed 
lines) and the green sample (thin dashed lines), while the right panel shows the finer color bins used to further divide the main red and blue 
samples. Bottom row: Mb versus redshift for red(left) and bluc(right) galaxies in the full DEEP2 sample. The magnitude and redshift 
ranges of the red and blue luminosity samples used here are shown with dashed lines, and properties of each sample are given in Table 2. 
The blue galaxy sample has been randomly diluted for clarity in this figure. 
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Fig. 2. — Coadded spectra of galaxies in the various color samples used here. The top left panels show average spectra for the main 
blue and red color samples, while the other panels show spectra for the finer color bins. The correspondence between restframe color and 
spectral type is clear. 
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Fig. 3. — Contours of constant correlation strength for the two-dimensional correlation function, i;{rp,n), for the main blue, red and 
green samples (top row) and various cross-correlation functions (bottom row). A 1 X I h^^ Mpc boxcar smoothing has been applied in the 
figure, though it is not used in calculations. Contours levels are 0.0 (dashed), 0.25, 0.5, 0.75, 1.0 (bold), 2.0 and 5.0. 




Fig. 4. — The projected correlation function, Wp{rp), for various color samples and cross-correlations. Errors are computed from the 
variance across the ten pointings on the sky. Corrections have been applied for observational biases due to slitmask effects (see text for 
details). 
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Fig. 5. — The clustering of the main red and blue galaxy samples, as measured by the u) estimator of lPadmanabhan. White, fc EisensteinI 
l|2007l ) (see §3.3 for details). Note that adjacent points are correlated. The solid line shows the best fit power law to these measurements, 
while the dashed line shows the result of a power-law fit to the Wp{rp) points in Fig. 4. There is good agreement between the two fits, 
given the measurement errors. 
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Fig. 6. — Deviations of Wp{rp) from a power law with rg = 4 h ^ Mpc and 7 = 1.8 for the main red and blue samples and the red-blue 
cross-correlation. The cross-correlation is seen to be closer to the blue auto-correlation on small scales, at rp < 1 Mpc. The geometric 
mean of the main red and blue Wp{rp) is shown as a dotted line for comparison with the measured cross-correlation function. 





Fig. 7. — Contours of constant correlation strength for the two-dimensional correlation function, ^{rp,iT), for the blue luminosity 
subsamples. A 1 X 1 Mpc boxcar smoothing has been applied in the figure, though it is not used in calculations. Contour levels are 
the same as in Figure [3] The strength of the fingers of god on small scales are seen to increase with luminosity in these samples. 




Fig. 8. — Contours of constant correlation strength for the two-dimensional correlation function, ^{vp, n), for the red luminosity subsam- 
ples. A 1 X 1 Mpc boxcar smoothing has been applied in the figure, though it is not used in calculations. Contour levels are the same 
as in Figure |3] Strong fingers of god are seen in each of the red luminosity samples. 




Fig. 9. — The projected correlation function, tUp(rp), for the blue (top ) and red (bottom) luminosity subsamples. There is little 
luminosity dependence seen in the fainter samples (left). Both blue and red galaxies show departures from a power law with a rise in the 
correlation function on small scales {rp < 0.3 Mpc) in the brighter samples with Mb < —20.5 (right). In the right panels the thin 
dashed line shows the Mg < —19.5 results from the left panel for comparison. 




Fig. 10. — Deviations of Wp{rp) from the best fit power law for the main blue and red samples for the brighter blue (left) and red (right) 
luminosity samples. A significant rise is seen on small scales for both blue and red galaxies in the brightest samples. 
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Fig. 11. — The clustering scale-length ro (left) and slope 7 (right) for DEEP2 galaxies as a function of median absolute Mb magnitude 
(top) and L/L* (bottom) for both red and blue galaxies. The values of ro and 7 for each luminosity sample are given in Table 2. Errors 
are estimated using the variance of power-law fits to jacknife sample s. In the lower panels we compar e the scale-length fo r red (left) and 
blue (right) galaxies found here at ^ ~ 1 with local results from 2dF jNorbere et al.ll2002bl) and SDSS llZehavi et al.ll2005l) . 
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Fig. 12. — The clustering scale-length, ro (left), and slope, 7 (right), of DEEP2 galaxies as a function of restframe color. The values of 
ro and 7 for each color sample are given in Table 2. Errors are estimated using the variance of power-law fits to jacknife samples. Dotted 
horizontal lines indicate the color range for each point; for clarity the color ranges for the main blue and red samples, shown as triangles, 
are omitted. The vertical dashed line indicates the approximate location of the minimum in the color bi-modality which separates red and 
blue galaxies. 
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Fig. 13. — The mean relative bias of each of the color samples compared to the main blue sample, averaged on scales rp = 1 — 5 h~ 
Mpc. Triangles show the main red and blue samples and diamonds show the finer color bins. Dotted lines show the relative overdensity 
as a function of color derived from the ^3 environment parameter of Cooocr et al. ( 2006), again normalized to the main blue sample color, 
which agrees very well with the relative bias measured from the correlation function. 




Fig. 14. — The quadrupole-to-monopole ratio, defined as Q{s) = ^2/(^0 — ?o)i as a function of scale for the main blue (left) and red 
(right) samples. Dotted lines show the expectations from linear theory for /3 = 0.65 for blue galaxies and = 0.52 for red galaxies, where 
/3 = fl^p /b and b is the linear galaxy bias. 




Fig. 15. — The two-dimensional correlation function, ^{rp,iT), for blue (left) and red (right) galaxies in the Millenium Run semi-analytic 
model mock galaxy catalogs ICroton et al.|[2006l) matched to the DEEP2 survey. A 1 X 1 Mpc boxcar smoothing has been applied in 
the figure, and contours levels are the same as in Figure[3] Almost no small scale redshift-space distortions are seen for blue galaxies, while 
red galaxies show strong fingers of god as well as coherent infall on large scales. 




Fig. 16. — The projected correlation function, Wp(rp), for blue (left) and red (right) galaxies in the Millenium Run semi-analytic model 
mock catalogs (solid line) and the DEEP2 galaxy data (dashed line). There is a strong discrepancy between the model and the data for 
blue galaxies on small scales [rp < 1 Mpc), where the model does not have enough blue galaxies. The red galaxy correlation function 
is too high on all scales in the model compared to the data. 




Fk;. 17. — The projected cross-correlation function for blue and red galaxies in the Millenium Run semi-analytic model mock catalogs 
(solid line) and the DEEP2 galaxy data (thin dashed line). The geometric mean of the blue and red auto-correlation functions in the 
semi-analytic model is shown for comparison (dot-dash line). On small scales {rp < 1 h^^ Mpc) the cross-correlation function in the model 
is significantly lower than both the geometric mean of the blue and red galaxies in the model and the cross-correlation function in the 
DEEP2 data. 



